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ABSTRACT 


The spatial distribution of power in a nuclear 
reactor tends to follow the neutron distribution because the 
energy resulting from nuclear fission is released largely 
at the point where neutrons induce fission. Leakage in an 
unreflected reactor causea both distributions to be hicher 
near the reactor center. The power at the center is limited 
by the maximum temperature which can be tolerated by core ma- 
terials with the result that outer regions operate at less 
than ultimate capacity. Ideally all regions should be opere 
ated at the maximum tolerable temperature, which would require 
uniform spatial distribution of power. Thus an important 
improvement in a reactor's power output can be achieved by 

aking the power distribution more nearly uniform. 


One simple, economical method of achieving uniform 
power consists of compensating for the decrease in power near 
the edge of the core by means of a reflector. If the return 
of neutrons from the reflector to the core can be adjusted to 
compensate for the leakage factor, the power distribution can 
be markedly improved. This method, applied to a reflected 
spherical homogeneous reactor, is studied in this papere 


The distribution of neutrons in a reactor is 
specified reasonably well by the diffusion type equations to 
be disevsced in Chapters I and TI. Multi-egroup methods of 
solution are used because or the importance of the dis- 
tribution of fast neutrons. Assuming two groups of neutrons 
was considered inadequate for the intermediate type of re= 
reactors which are studied. Multi-group calctlations are so 
time-consuming that the use of computing machines is essential. 
Chapters I and II develop techniques for use with M. I. T.'s 
digital computer Whirlwind I. 





ABSTRACT 


(continued ) 


To simplify the problem, a basic reactor is 
selected witn\highnlyvenriched fuel, sodium coolant, beryl- 
liun moderator and reflector. The reactor is assumed to be 
homogeneous, and: inthe steady state. The effect of reactor 
poisoning is not considered. Spherical geometry is assumed. 


It was found that by suitable adjustment of the 
naterial proportions in both core and reflector the average 
power could be raised to within 16% of the maximum power. 
Results are presented and discussed in Chapter III. 


Thesis Supervisors: Thomas H. Pigford 
Melville Clark, Jr. 
Titles: Associate Professor of Nuclear Engineering 
Assistant Professor of Nuclear Engineering 
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ieutron leakage from an unreflected reactor results 
in a lower concentration of neutrons toward the reactor edge. 
R3y virtue of the fact that the power resulting, from nuclear 
fission induced by neutrons is largely liberated at the 
point where the fission occurs, the spatial distribution of 
power also is depressed near tne reactor edge. This imposes 
@ limitation on a reactor's power output. The maximum power 
censity at any point is limited by temperature considerations. 
Outer regions, With lower power density, operate at tem- 
peratures lower than tneoir ultimate capacity. Elementary 
considerations indicate that the average power in an 
unreilectec spnerical reactor is only one third or the 
naximum power produced at the reactor center. Tous, if a 
unitorm power distribution could be achieved, an important 
inercase in reactor output could be realized. 

opatial variations of the proportions of tne mate- 
rials in the core, together with use of a reflector, is one 
Met20d OF isproving power distribution. Another is the use 
Or &@ poison to depress the neutron concentration near the 
reactor center. These methods have the disadvantages of 
encreased complexity and cost. Aa simple, economical alter- 
native, the use of neutrons reflected back into the core 
from the reflector to compensate for tne vower depression in 


the outer core, is investigated here. 





The diffusion equations to be discussed in 
Chapters I and II specify reasonably well the distribution 
of neutrons ina reactor. Since neutron energy varies 
widely and since the coefficients of the terms in the dif- 
fusion equation denend on neutron energy, some simplifi- 
cation 1s necessary. One method, applied to thermal re~ 
actors, assumes that neutrons are slowed to thermal energies 
immediately after birtn. The coefficients of the terms in 
the diffusion equation are then determined at thermal 
energies, and a gimple analytical solution is possible. 

However, in this probiem, the number of neutrons 
with energies greater than thermal must be increased in order 
to sunply the reflector witn fast neutrons for moderation and 
reflection back into the core to compensate for leakage. 
Therefore, the distribution of higher energy neutrons must 
be known, a&nd one=rroup calculations are inadequate. 

Considering two groups of neutrons—-—fast and 
Slow--would give some representation of fast neutrons which 
leak from the core to the reflector, where they may be 
moderated and returned to the ccre as slow neutrons. Howe 
ever, better solutions are possible if more than two groups 
are considered. 

Multiegroup calculations require solutions of dif- 
fusion equations for each energy group of neutrons con- 
Sidered and are very laborious and time-consuming. The use 
of a high speed digital computer serves to minimize both of 


these problems. In Chanters I and II, diffusion equations for 





four and six energy groups respectively are broken down into 
a set of simple calculations which a comvuter can perform to 
Serve the equations. The actual programmed set of ine 
structions is used with M. I. T.'s digital computer, 
Whirlwind I. The basic bs ibe is applicaole to other 
digital computers by a suitable transformation of the 
snecific set of instructions. Whiriwind I requires less 
tnan a minute to solve the equations, but a period of from 
two to three minutes is required for the comnuter to transe 
form the input of decimal numbers and interpreted instruc- 
tions to its own binary system. This conversion time can be 
reduced by rewriting the program in the basic Whirlwind code. 
The Whirlwind I computer will perform this function auto- 
matically and will cut a new tape containing the converted 
programe 

In this paper a basic reactor composition and cone 
figuration are chosen, and the proportions of the various ma- 
terials of the core and reflector varied in order to dee 
termine their effect on the power distribution. The funda- 
mental reactor type studied was cne considered suitable for 
high output mobile reactor application. 

The requirements of this type of reactor led to 
the use of core materials capable of withstanding high 
temperatures. <A liquid metal coolant was chosen in order to 
avoid having to pressurize the core. Of the liquid metals 


available, sodium was chosen because of its superior heat 


2 Cask 





removal characteristics and because of its comvatibility with 
other materials at high temperature with resnect to corrosion. 

Once sodium was chosen as the coolant, stainless 
steel became the best material for cladding in view of its 
high temverature strength and resistance to corrosion by 
sodiume 

A dispersal fuel element of uranium dioxide ina 
matrix of stainless steel was chosen to meet high temperature 
requirements and to minimize radiation damage. Galculations 
for a fuel element of this type have been me 

Beryllium was chosen as a moderator and coolant 
because of its high moderating ratio and good thermal pvroper- 
ties. 

The average operating temperature of the reactor 
was taken as 1000 °F, which corresponds to a thermal neutron 
energy of 0.07 ev. 

For mathematical simplification, spherical geometry 
was assumed. The reflector thickness wag uniform. Both core 
and reflector were assumed to be homogeneous and in the 
steady state. The effect cf reactor poisons was not con= 
Sidered. Thus an actual reactor would require additional re- 


activity to overcome the effects of fission product poisons. 





“Proceedings of the International Conference on the Peaceful 
Uses of Atomic Energy," Volume 9, page 561, C. Ee Weber and 
He H. Hireeh. 
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CHAPTER I 
Method of Finite Differences 


Le N. Hoover 


= 


This chapter will develop a program for solving 
multigroup neutron diffusion problems in spherical geometry 
on the Whirlwind I (WWI) digital computer using a finite 
difference approximation to the age diffusion equation. The 
solution will give the power distribution, the neutron flux 
distribution in each of several energy groups, and the 
critical radius for a reflected spherical nuclear reactor 
having a homogeneous core and a homogeneous reflector. The 
program will be developed for neutrons in four energy groups. 
By suitable modification of the basic program, the number of 
energy groups and the number of regions considered can be 
increased. Four energy groups represent the maximum number 
of groups which this program will conveniently handle using 
only the internal storage of WWI. To increase the number of 
energy groups considered will require using auxiliary 
storage, which has the disadvantage of increasing the 
machine time required for a solution. The four~group 
program requires approximately three minutes of computer 
time for a complete solution. For a parameter study, the use 
ope a larger number of energy groups was not considered to: be 
ee since very little improvement in results would 


be achieved by the use of more groups. 
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I. sge Diffusion Equation 

The age diffusion equation provides a good de- 
scription of the neutron behavior of a large class of 
nonthermal reactors. The following development is concermed 
With reactor assemblies that are large compared with the 
neutron mean free path, and in which the neutron energy is 
reduced primarily by elastic collisions. The method accounts 
for both fast and thermal fission, for tne fission spectrun, 
and for inelastic scattering. 

The basic equation used is the age-diffusion, which 
gives good results for the case of weak capture anda 


2 
cross section that varies slowly with energy, 


page) + 2, (u)f(u) = Seba + S(r,u) (ee) 


Tne mies - Satu) represents the excess neutrons entering the 
lethargy interval over those leaving, assuming, that the slow- 
ing down density, q, can be treated as a continuous function 
of the lethargy, us» This excess number of neutrons is 
balanced in tne steady state by tne leakage of neutrons out of 
the region on interest, which is =p7V@(u), where Z(u) is the 
flux per unit lethargy interval, and by the neutron adsorption, 
2, (u) G(u). The term S(r,u) represents neutrons entering 
lethargy u ait a space point r as a result of production by 
fissions, plus inelastic scattering from a higher energy. 
There are advantages to expressing Equation (1.1) 


in terms of the slowing down density, q, rather than the 








=$S 


flux, @. The relationship between the flux per unit lethargy 
interval and the slowing down density in the case of slowly 
varying capture cross section and weak absorption Lk’ 

Z(u) = Eater = (1%) 
For energies below the sowweeenergy, the same number of 
neutrons must sicw down past each energy, except for the ef 
fects of leakage and absorption. Hence, for reactors with 
§$mail absorption in the slowing down range, the slowing down 
density tends to be a less sensitive function of energy than 
the flux, %, which will vary more rapidly than q as P* yeas) 
variese 


In terms of q, Equation (1.1) becomes 


-0v/ eats |+ Spr 2. Suu + 9(zem) 


In Equation (1.3), uw is a continuous variable. 
For multigroup calculations, Equation (1.3) 4s integrated 
over the lethargy range from u, to uo, corresponding to the 
width of the i'th energy group. 2 ,{u) is considered as 
inderendent of position and brought in front of the 
Laplacian cperator, Vv ~ The agséeediffusion equation for 
lethargy group i then OY a 


D uw) OUT 


( D. (2a q) = 
UiCFE,Y 4 re “4G Ms ave “4 


og q, 7 + fotzyuran, (1.4) 
U; 





Lie 


Ue 

where Us =f du. The averages are With respect to u in 

u 
Al Nw 

tne Groupe cre is the degradation of neutrons into the 


i'th group from the grouv of next higher energy, and qe is 
the desradation out of the i'th group. The approximation of 
replacing the averages of products apvearing, in Equation(1.4) 


by the product of averages yields 


2 
r up) er (222) 
> q). a8 : (q) 
F2s/, ms i ave Gz 5 i ave 1 ave 
IN | 
rc ec ee) S(a,u)du (1.5) 
Ui. 


Mate precedure is justified if either gq or the cross sections 
Go not change radically over the group. 

To solve Equation (1.5) when ae and 5S are ixnown, 
it is necessary to postulate an additional relationsnip 
between (q)4 ave: alt and ee The simplest method would 
be to assume a linear relationship for q across the lethargy 
Group and to take an arithmetic averare. However, with 
only four lethargy groups, the large range of lethargy per 
group tends to make the linear relationship inaccurate. A 
Suitable method for use with a small number of energy froups 
is to use the relationship 

(a)4 ave = Sane (1.6) 

where f is a number chosen by trial and error for each 
group and region of the assembly so that f times tne spatial 
average of ae gives a correct value ror the spatial 
average of (q); aves The solution is not sensitive to the 


choices of f, and it is usually possible, on the basis of 





Hie. 


experience, to make adequate choices the first time. From 


the problem solution, the values of (q) are plotted 


i avg 
versus léthargy as in Figurele. A smooth curve is faired 
through these values of (a)4 aye to represent the actual 
process of continuous slowing downe From this curve, 

refined values of f can be determined and used to improve the 


solution, if the initial estimate was greatly in errore 


II. Average Nuclear Properties of Material in 2 Lethargy 


Group 
Equation (1.5) expresses the age diffusion equation 


for lethargy group 1 in terms of material properties averaged 
over the ifth groupe Within each group, the neutrons are 
assumed to diffuse without energy loss until they have suf- 
fered the average number of collisions which would be 
required to deerease their energy to that of the next lower 
groupe At this point it is supposed that the neutrons are 
suddenly transferred to the latter group. This process is 
supposed to continue while the neutron energy is degraded 
from the group of birth energy to that of lowest (thermal) 
energye 

Since the energy spectrum of the neutrons ina ree 
actor covers a continuous range from thermal energies to 
about 10 Mev, the croupediffusion concept is evidently an 
approximation. This is partially overcome by assigning ap-= 
propriate average values for each group to the various 


properties, such ag the cross section and the diffusion 





oe 


C@mliiaci@n. of tae cor@and the refiector. 

Omit: USO Cee 1 icrent. Thevwavergge diffusion 
CO@ii 1cheit, (04) {or Gncway group 1 is derived from a con-= 
sideration of the fust neutron current density. The 


mL Ini §2en O01 De is 


ail LD . 
file) 7 B(x, 2) aE 
iD = =o Ae (250) 


sk Ey 
Jn, V Meee 


It can be seen from this enuation that Ds has the same value 
aot Giiterent spatial points witiuin the region being con- 
sidered only if @(r,E) is separable into a function of r 
multiplied by a function of FE, i.e8:, if the variation in the 
neutron flux spectrum from point to point in tne medium can 
be neglected. This assumption is usually made for Simplicity. 


Hence, with D(H) replaced by ~—_--___, the expression for Ds 


PBs GS) 


becones i 1 
(Fe ee 
is 
f. 1 ge) aE 


Eo 


Kor the asymptotic case where EX<&E,, which is assumed to 
apoly +or all energy Groups considered here using Be as 


the moderator, and for the case of weak absorption 


Be 03) - af) 
~~ EE Z,(E) cas €Z,(u) 


In lethargy notation we find that Equation (2.2) becomes 





Bs 


e 


pa L sretthtier 
| sake 


u, 


eM! a, 
= _ OF OF Ztr2s)s ave = (3 5 Z¢rZ)4 ave ‘1 ave 
a F28)4 ave (FZ5), ave 44 ave 


The last step involves the approximation of replacines the 
averares of products by the product of averages. This is 
justified if either q or the cross sections do not change 
radically over the groupe Hence, this is a criterion which 
must be considered in choosing the energv subdivisions of 
the group. Thus the average diffusion coeffisient for 


rroup 1 is given by 


a 
p, = Lo f4sztr)s ave (2.3) 





be Total Flux Per Energy Groupe fhe total flux 





in energy group 1, NV,, 18 defined as 
u 


a 
NV, = ) Badan = fern 
u U: f2,) 


Use of the same approximations as for the diffusion coef- 


a! 
ficient gives the result, whereU, = f au 


1) . 


ia | 
NV, = Us x5 4 ave (aq), ave (2.4) 





- 8. 


Or 


NV, 
(Qs ave = -—_—y (24a) 


MES q)4 ave 


c. Average Absorption Cross Section Per Group. 





The average absorption cross section for energy group i is 


defined as 
ae = (u)q(u)au 
ee JZ (u)g(ujau - 2, §e.tay 


EOE ———— § 
42 u :}du 
J B(u)du y! : AAS 


When the same approximations as before are made, the average 
absorption cross section for energy group i becomes 


Za, 7 (2 /55)4 ave 


y G2Za°1 avg (205) 


III. Spatial Intezration of the Group Diffusion Equations. 





Equation (1.5) describes the neutron behavior in 
energy group 1. Substitution in (1.5) of (2.5) plus the 


following relationshivs develoned previously 


i are (1.6) 
: (F638 try ave 


ID ae ea ec 
1 -~ ( 1. (2.3) 


) 
CE da)4 ave 
NV 
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(2.4a) 
= as a ae 
U ee 

1( F254 ave 


(q)4 ave = 


gives 





=~ Gow ( 


Dy ONY, i“ (ay fi, {2 | NV 


CRS 8)4 ave) 


J 
ce a, + ans ical (3 ils) 


For spherical geometry, the following substitutions will 


simplify (3.1) by eliminating the Laplacian opgrator vee 


Ly - rv, (3.2a) | 
1 Our OUT 

Q ~ Sia (3.2b) 

i3,, = rn J5trpat)au (3.2c) 


A new notation is har Let i refer to the ifth 

energy group as previously and the subscript n refer to the 
n'th space point located at a distance r,, from the origin of 
the spherical coordinates. As before, the source term th 
representg neutrons entering group 1 at radius ry as @ 
result of inelastic seattering and fissions in all grouns.e 


Use (3.2) in (3.1) to get 





2 ip i 4 
= Ds d 5 n+ Pay re 1 = +P, va Qn OUT 4 on 
ar ae a 
Pata i i (3.3) 


The space points are chosen to divide the spatial region 
{nto a set of subintervals of width Or. The finite dif- 


ference approximation . . 


gi ip ie — 2 ip 
{a pda i > ntat Fy = 2 : (3.4) 
(Ar) 


and Equation (3.3) gives the recurrence relation 








he & 4 di a 
Put. = Ky ms i oe - “I, (3.5) 
where 
LiL - Dy 
b= Ja, + 70, at (3.5b) 
§28'4 ave 
2 = ro | 
tr, = eae [stg 77 + 4s, | (3.5¢) 


Equation (3.5) can be verified by direct substitution of 
Equation (3.4) into Equation (3.3). The above substitutions 
follow a method developed by Ehrlich and iat, 

ae Boundary Conditionse The boundary conditions 


for the numerical solution of Equation (3.5) are that 


ips 0 at r= 0 and also at the extrapolated end point out- 
Side the assembly. Another condition is that the neutron 
flux and current density must be made continuous at 
boundaries between regions of different material properties. 
Suitable choices of Zr must be made for each region so that 
the boundaries will occur at space points. At a boundary 
point, M, located at radius r- ry between two different 
media, these boundary conditions give (primes refer to the 


outer renton) 


|» 


i ‘FH, hor cantinurty orf ius, 


Jus Jy, for continuity of current density, 


where 





: 
7: 





_ 
_ 














: OO 





. _ as OO 


7 : -_ 
7 ow - : 
> A A ee 





gop 20M) - p(eF +r 9) 
= ay ~ ais ar 
It can be shown that (seé Appendix C) 


x ¢ 
(aha | ee (Pusar ~ Fu - Ar) - = Py (3.54) 


oa (igi Ar T Fy . Ar! 
Therefore, the boundary condition for neutren current 


density, J= J', is 


a: rl * 
3 Py ' 6 (Pu¢ar Fy -pr? it (Fy iar 7. By ~ As? 


« > 2 : 
-- | $ Pit ZF yt Aget Fly Apt) - ty 5 


Hy re? one (3.6) 


b. Method of Solutione The recurrence relation 








for 1p Equation (3.5), will be solved using a procedure 

suggested by Ehrlich and Hillel te" This procedure is 

particularly convenient for use with the WWI digital computer. 
To solve Equation (3.5) subject to the above 

boundary conditions, new varfables Cc, and Be. will be intro- 

duced. These variables are related to the homogeneous 

(In = 9) solution of (3.5), A,» by the defining relation 


Ships 
n 
1% 


- Anti. ana S A 
~ntics= An Pa < ye 


When the.above expression for & is substituted into (3.5), 





with be - 0, it is found that 


o¢ 


aa 


£22 = on (3.7) 





as oe 


The recurrence relationship for tne(@'s is 


2 = ee I (3.8) 


WilCn Con We™weriivea by direct substitution of the définition 


of @ into (Ceomremecomoination of Equations (3.5), (3.708 


and (43.6) gives the followins relationship 


a ney + P., 
iy OS (3.9) 
Kel 
This is proved as follows: 
Substitute (3.7) and (3.8) for &and into (3.9) 


and combine terms. This gives 


-. = > 
eG F, + Baal Fak a nels I, 


By Equation (3.5), the right side becomes F,_j. 
This is then Equation (3.9), when n is replaced by n¢l. 

Equation (3.7) is used to solve for the &'s, 
Startvlive With Xo =k, since A, = O and therefore,from the 


definition of &, &, =O. The next step is to use (3.8) to 


solve for tie 3 "Ss, Starting with 2, = £2, ranally (4.9) 


it 
is solved for tne “'s from tne outside to the center of the 
spmerical reector ascenbly, starting with F,, - O at the 
extravolated end OUN b. 

The ecuations for & and @ have the same mathe- 
macurcal 1Orm DUG Giiferent constants in regions of difiterntwase 
material properties. at the boundary between regions, « and 


QO are Single valued in order to maintain continuit: of flux. 


Thus Equation (3.9) is unchanged at a boundary. At the 





ie 


boundary point wi, et radius ro, Separating two different 
media, the boundery condition for continuity of current 
density can be satisfied by using Equation (3.6). The sub- 
stitution of the recurrence relationships for & and © into 
(3.6) leads to the following relationships across boundary 


point i (primes refer to the outer media): 














t wo ial Lue! 1 
on sl : sfir Xp (5-10) 
‘ - es Ar' ‘ 
GB i pl, + (2 1g “Oe ‘ M 
oe 
a AGE el (3,198 
sAr &, 
where 
ee es 
D 
jmeeae’ = LN! 22 Ay! . 
a) — ene —_ — PRB 4%, le 
: 2s Ar 20 y 4 SR (3 
Pee Las ior 1 eae 
xk = pk + fl + BP, 7s ( - Br, 


Cc. liverral Group. Thettbermal neutrons are ase 
sumea to be in a group of essentially constant Wier. It 
is assuuea that neutrons enter the thermal group only by 
elastic Scattering Prom the group of nOxt higner energy. 
the tnerical group is treated the saine as other groups except 
fOr chen@ins tHe definition of 1 Vin (where tne subscript T 
Stends for the thermal group). ior the thermal group 


*% for the case considered here with 15 space Doints in 
the core so that Ar= R/15, where R is the core radius. 
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NVn g Lr 

de Energy Groupse Energy Group I is chosen from 
10 Mev to 0.1 Mev (u= 0 to u & 4.606). Glasstone- plots 
the fission neutron energy spectrum, and it can be seen from 
this plot that very little error is introduced by assuming 
that all neutrons released in fission are born into Group I. 
Groups Il and III are chosen to give approximately equal 


lethargy widths down to a thermal energy corresponding to 


1000°F, 
The energy groups will be as follows? 
Group Lethargy; u_ Energy, E Lethargy Width 
a _ 0 w 4,606 10-0,1 Mev 14,606 
ie 4.606 = 10.81 1006-9 .2 Kev 64204 
III 10.81 ~ 18.80 200~0.070 ev . tages 
r 18.80 0.070 ey (1000°F) 0 


@e Inelastic Scatteringe Provision is made for 


considering inelastic scattering into an energy group by the 
definition of the term | S(rgu)du in the age diffusion 
equation. This term ft absiinnt to be composed of fission 
neutrons born into group i plus those inelastically 
scattered into group i from higher energies. Since a 
neutron must possess a minimum of about 0.1 Mev of energy if 
it is to be involved in an inelastic scattering process ,° 
inelastic scattering occurs only in Group I. The group 
widths are large so there will be no inelastic scattering 


across groups. This can be seen from plots of the energy 
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distribution of inelastically scattered neutrons as given in 
the Reactor Handbook, Reference (7). 
IV. Programming for the Whirlwind I Comnutere 

Whirlwind I (WWI) 1s a digital computer capable 
of performing simple mathematical operations (addition, SU Dee 
traction, multiplication, and @ivision) with extreme 
rapidity under the control of a set of instructions, called 
the program, which is stored in the computer. A detailed 
account of WWI computer operation is given in Chapter II, 
Section II, of this thesis. For additional information, 
the reader is referred to Reference (1), which completely 
covers WWI operating techniques. 

a. Numerical Onerationse Section III of this 
chapter develops the relationships which are used in 
programming the WWI solution. The steps of the solution are 
outlined below: 

(1) Start by assuming the spatial distribution 
of fission neutrons born into Group I as a result of fissions 
in all groups. (This assumed distribution will be improved 
to the correct value by successive iterations.) 

| (2) From the assumed spatial distribution of fige 
gion neutrons into Group I, which is J sewn, use (3.2c) 


to get es. Ur 


oe = =p JsCrewau, 


Ur 
(3) Oaleulate the following constants for Group I, 
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Mein oouaeions (5.5) and (3.12): 


= i 
warn ae 
PiU; i 
§2s 1 ave 





2 oe 
‘ D 
p= Ca Aes 
esitAr 3r 
1 (GaN 
k 
Di 
’ ' t 
ary, ff © ) oe 
(4) Calculate co starting, with ‘Xn Ik and using 
urd it FS As \ — | _ i - ra + 
PemmelOn (4.47). Xe] K a The value of k must be 


correct ior the resion considered. at the noundary betweén 


6Oro 2°C werl@®ctor, where’n = 15, boundary condition (3.10) 





: es 
is used a a ee 
’ Kx 16 ig eye a ; 

(5) Calculate Ig, Starting with “8,2 “I, and 

Mot OGe Ge tLOn (| >.C) , B= —— + I, where 
EN 
a ae cS ts forl#n‘515 
4H D Yh 
I 
qi = 0 for 16 =n # 25 


me o@e boundary, n = 15, Equation (3.11) apvlies, 


AN 1! Aeon Br», 
aS sfAr M15 





E pel 
Bo. > ORG as (2 + 6ric 


(6) Solve for 1s starting with ‘Yo, = 0 and 


WOr.cine tow@rd the center of the spnerical reactor assembly 





a>... 


using Equation (3.9), 
ty 2 TF, - as o 
Snel 
\ Il, 
(7) Calculate the inelastic scattering, rage 


from Group I into Group II by use of the relationship 
= ii 
i Me 


n 

where X is the probability that a neutron scattered in 
elastically in Group I will enter Group II. 

(8) Repeat successively steps (2) through (6) for 
Groups II, III, and T. Note that there are two changes to 
the procedure used with Group I: First, the input term, 
te is composed only of inelastically scattered neutrons 
entering Group II (since all fission neutrons are born in 


Group I), and a ig zero for Groups III and T; second, the 


calculation of +I, by Equation (3.5e), 
iy _ (Sn)? 5 +aeg 
n= Dy Q n 


1—1 
ineludes elastic scattering, q.. aa out of the group of 


next higher enerzy.e 
(9) When “A. has been calculated for all energy 
I 
groups, calculate a new source term, S,, of fission neutrons 


for Group I from the relationship 
: 4 
eo z (Y2e)4 ave Pn? 


where (vy 2), ave is the average value in group i. fhus, 
Is. is proportional to the fission power distribution in the 


reactore 





=} O= 
(10) The new input tern, “o is normalized and 


used as the basis for a new iteration, which repeats the 
above procedure starting with step (3). 

bd. Iteration Procedure. The itcration procedure 
has been divided into two steps: (1) relaxing from the 


assumed input of fission neutrons, Is to the correct 


Nn? 


spatial distribution; and (2) changing the radius, R, of 
tie NeW Meactor core to obtain criticality. 

Let the superscript m indicate the number of the 
iteration, with m = O indicating initially assumed values. 


Let a circuvflex indicate values which have been normalized. 


I 


Step (1) corrects the spatial distribution of “S, 


2S tollows:; 


All iterated input terms ‘s,” 


froin the assumed “ean 


amco Group I, Stara 
are normalized with respect to an 

; ; : ie meal . 
arbitrarily chosen core spatial point, N, where Sy 4 is 


normalized to unity. A normalizing factor & is defined such 


evele 
© mel 
Sy ah 
Se “a 
S71 Ory 
and AN 
, ig m yls m 


wnere WF = the chosen spatial point for normalization to unity 
m= tne nunber of tne iteration 
5." = the calculated input term to Group I et spatial 


jayenhavigee! 





Ke 


IA mo 
Sh - normalized source team 
IA m=-l - ee 
Sy ~ 1 = source term of previous iteration normalized 


to unity at space point N 
To determine when to stop iterating, the normalized 
LA 
input term, ar is compared at each spatial point, n, with 


the preceding iteration (see figure below). It is 




















A tl= 
| se 3 Is 
ae: Core Power wy 
ae Distribution = 
aa —— a eo 
“Sh at SS ca 
normalization | 
DemneL | ! 
| | 
O it a DD so we seteeeeninre eee es 


specified that at each spatial point the ratio 38 

must differ from unity by no more than a small value, € , 

to permit the prorram to proceed to the next step. Whenever 
the value € is exceeded at any space point, the program “ill 
reiterate. This is equivalent to requiring that the solution 
for the core power distribution must have converged suffi- 
ciently to permit the solution from one iteration to lie en-= 
tirely within an envelope determined by around the so- 
lution of the preceding iteration. 

The criterion for choosing between iterating or 
continuing to the next step is based upon the ability of #wl 
computer to choose one course of action for a positive 
number and a different one for a negative number (conditional 


control transfer). Provision is made for the criterion to 
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operate, irrespective of whether 5, “/S, is greater or 


less than unity, by the following scheme: 


G tl 
(a) When Sy is greater than unity: 














m 
Pn 
g ml 
l= +€ <0, not acceptable (reiterate) 
~~ i 
Sy 
“N 
g mel 
1S lt : +¢€ 70, acceptable (proceed to the 
on next step) 
3 me ] 
(b) When — ig less than unity: 
~ 
Sy 
g_ m1 
n ~l1 +€ € 0, not acceptable (reiterate) 
“ 
Sy 
5 mm]. 
n ~l +€ ? 0», acceptable (proceed) 
A ih 
Sy 


Step (2) of the iteration procedure involves changing the 
core radius to get criticality. The procedure is as follows: 
Step (1) above has required that a plot of the 

input term s must reach a stable shape before proceding to 
step (2). In step (2) the last calculated value of y by 
(Squation 4.1) is used to test for over or under criticalitye 
A & less than unity means that the number of fission 
neutrons entering energy Group I increases from iteration to 


iteration and, therefore, that the multiplication factor, k, 
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is greater than unity (and vice versa). For the critical 
> shape stabilized 
condition, k =~ 1, the,source term remains constant for suce 
cessive iterations, and do therefore is unitye 
If J is'not within a small value € of unity, the 
program goes through another iteration using a new value of 
the core radius, R, chosen to come closer to criticality. 


The new radius for iteration mis 


Re = Y mel oe eat 


Thus, the core radius is increased 1f k <1, and vice versa. 
The criteria for determining whether the radius is correct 
or not is: 
(a) When J is ereatver than unitv; 
ewe 1 ao not acceptable (modify radius 
and reiterate) 
Em Xt 1 70 acceptable (proceed to the 
next step) 
(b) When % 48 less than unity: 
y- 14 <0 not acceptable 
ySIt+e 76 acceptable 
When the iteration criteria of both step (1) and 
» twp (2) are satisfied, the program proceeds to the final 
step of calculating and printing the core power distribution, 
the critical radius, and the neutron flux distribution in 
each energy grouvde 


Ve. Flow Sheet. 


The program for the solution of the reactor problem 





has been described in Section IV. A flow sheet will be 
developed to show schematically the operations discussed 
in Section IV. The following conventions will be used? 
A square box standing on one 
edge denotes a sgubroutine for 
calculating the quantity within 


the boxXe 


Circles enclosing numbers are 
used to specify corresvonding 
parts of the flow sheetse The 
(6) numbers are called “addresses,." 
The circles are called connectors 
and are put at those places ina 


flow diagram to which reference 


will be made at some other pointe 


A rectangular box denotes that the 
— quantity inside is to be calcu» 

lated or the operations within 

are to be performed. Frequently 

the formula for calculating the 


Quantity is given. 


This symbol indicates a conditional 


transfer to a if A> B and to b 


i fo —G.. 
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Subscripts and superscrints will have the followe 
ing definitions; 

i denotes an energy groupe 

n denotes a space point at radius rye 

m denotes the number of an iteration, where m= & 
denotes values initially assumede 

N denotes an arbitrary space point, n= N, chosen 
as a basis for normalization procedures. 

> the symbol a+b means that the quantity a replaces 


the quantity b in the storage register in which 
b Is stored. 





22). 


FLOW SHEET rOR FOUR GROUPS 


Main Program 
















PoteGr.. 1 
constants into 
operating, registers 






Stored 
Constants 
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oc LON DEEL POR ec Ouk aiCues 2 


Group Constants Subroutine 








boundary 
cond2bion 





=ou = 
FLOW SHEET rOR FOUR GROUPS 


Z 
GZ. Subroutine 
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PLO. Meet FOR 2OUR GROUSS 


B puoroutine 











] Ar'hs dr’ eee 





ree Tee form of 1, comes from combination of 
fee ec), ond (425e) with Ar R/15, 
radius. (There are 15 core evace voints.) 


Boquations (iso 


where Rn = core 
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FLOW SGEET FOR FOUR GROUPS 


a Subroutine 


i 
1£n#25 













in 16 | oe - Fy 
= 2 “ R + (n-15) (4r_) 











Memorv Assignment (See WWI program following) 
Flad ” Contents 
pl storage of numbers used in numerical operations 
m2 Temporary storavre 
bl R, (YV2p)qz, VEe)riz, (YZe)n, Tr (4r*) 
eee O64,,)', € ; (v2,.), 
11 Used to indicate number of iterations 
ol Group I constants: D, Dt, fr, r°, Ge 5 ee 
(E20) 4 ave? Za? Ze!) L/y, 
Se Group II canstantsa: (Same) 
C3 Group III constants: (Same) 


*" A pLAD is @ floating address used to locate a particular 
set of storage registers 








if. 
f2 
t3 
fh 
Zi 


Ze 


23 
att 
x1 
ll 
re 
ae 
r4 
r5 
r6 


r7 
rg 


r9, rlod 
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Contents 
Group T constants: (Same) 
TJ 
Assumed source Sy 
Calculated source term from fissions in all 


T 
1 - 4m 
groups Zi Uw oF) - Sy 


e] 


a 12n #25 
ce Ll4n25 
a 

PF, 1lSn225 
st 


Operating register for storage of registers cl 
when group iis being used 


a 
Temporary storage for <x 


Temporary storage for F 


Temporary storare for Fn 
If 
Temporary storage for inelastic scatt. (3 an 


Group constants subroutine 


a subroutine 
12 subroutine 
n 


Cp subreutine 


Ip. calculat&ion stored in fl 


I 
Fn calculation stored in f2 
alee 
Bn 


ik 


calculation stored in f3 
F, calculation stored in r4 


Calculates (1s jan inelastic scattering and 


stores in xl 








r13, 810 


yl, sll 


X2 


al and yl 
az and y2 
a3 and y3 
at and y4 
38 


89 
s10 
sil 
912 
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Contents 

19 subroutine 

Iterates for correct spatial distribution of 
flux 

Changes radius and reiterates until criticality 
is achieved (both r13 and ri4 will verform a 
maximum of 50 iterations before printing out 
the answer 86 as to avoid a waste of computer 
time should the solution be diverrent) 
Caleulates “30 from the assumed 5, where 

Is’ i nAis.) 

Control for Group I calculations 

Control for Group II calculetions 

Control for Group III calculations 

Control for Group Tgealcutations 


Calculates X= 1 where N is spece point 10 


Syn 
Normalize 3." = wees 
(See r13) 
(See rl4) 


a uf 
Puts Bn in ¥55 ro, ©7, ro in place of Bne 
Prints core power distribution, flux distri- 
bubwvons, Critical radius, oa and number of 


iterations 





fe TAPE 374-223-008 
46) 


5 
t=epe 
pi,+1.0 
+720 
+13.0 
+50.0 
+0.0 
+20.0 
11,+0.0 
+*1.039 
+0,000196 
40.522 
+0.0 
41.05 
40.12484 
+0.,03734 
$1,+1.33 
+1.1 
+0.9 
$2,+0.0 
P1,+0.0 
£3,+0.0 
Z1,+0.0 
23 9%0.0 
x1,+0.0 
fav21+4 
itsz14+20 
iLdvbi 
imrzi4+i8 
imrzi+10 
ladpi+te 
imrzi+28 
icspi 
idvbi 
icazi+24 
fadzi+c4 
idvbi 
isco 
1tsZ24+30+¢c 
fcabi 


imrsi 
ROnnennaa 
ie BEC 
1@c01r 
imrbi+10 
1t82Z3+20 
1%s23+30+¢c 
icabi 


$caZzit36 
mrt 


+2.0 
+5,0 
$14.0 
D2,+0.0 
b1,+60.0 
+2.0 
+0,.0 
+0.621 


41.528 
+0.16119 
+0.03813 


c4 ,41.190 


40.1282 
+1.31 
+1.10 
+1.44 
DUO 

Li gee 
TO 
DITTO 
DITTO 
DITTO 
jmrzi¢8 
icapi+20 
itszi+oo 
idvazi 
iadzi+14 
itsz1i428 
idvpi+e 
tdvzi+eo 
tadet 
itsze 


1TSZl4O+-e 


idvzi+20 
icrg 


idvp1+28 
itees 
icaz3+ 
iceazi+cec 
{dvbi 
isec5 
4dvpi+28 


imr21+40 
itsz23 


wee 


WHIRLWIND PROGRAM 


A HOOVER | 
+3.0 4.0 +5.0 +6.0 | 
+9.0 +10.0 +141.0 +12.0 
+15.0 +16.0 +0.0 450.0 
+0.0 +0.0 +0.0 +0.0 
+0.001073 +0.01444 +40.052472 
+0.0 +0.00 +0.03 +0.000219 
+0.0 +O.0 c1,+2.209 | 
+0.68 +40,02158 +0.06884 
+0.0 40.2171 ©2,40.5756 
+1,08 +0.06535 +40.13231 +40.000523 
03 ,41.104 +0.537 +1,66 
+0.12968 +40.008951 4+0.000522 
+0.544 +1.0 +1.0 
+0.0354  +40.003914 +0.0 
41.29 +1.26 +4, 22 +1.20 
#1.05 a ove $0.95 40.92 
$1.45 Be eae 
28ri +0.0 
Sr) +0.0 £2,40.0 DITTO 4Sri +0.0 
48 +0.0 £4,4+0.0 DITTO 48x) +0.0 
HO, 40.0 22,40.0 DITTO 46x] +0.0 
48x, 40.0 24,40.0 DITTO 48x) +0.0 
48ri 40.0 vi1,itad54r 10az1i416 
fadzi412 itszi418 icazi¢+2 idvazi 
imrpi+2 idvpi+4 imrbi+10 idvzi+20 
Lea 01 idvpi+28 imrb1 Ldvpi+23 
fjadpi+2 itszi¢24 icagi+16 idvzi+6 
4tszi426 imrb1410 imrbi+10 idvzi+e2 
icabi+10 idvpi¢+4+  idvbdt iadpi 
ite leazgi+22 imrzi+24 itset 
jadp1 imrpi+2 idvpi+4+  4mrpi4+10 
Lladt itszi4+30 isp0 P2,1¢a25r 
isc4 ier13 iesp1 Ldvz2+ec 
ict5r lespl imrp1+25 
ddvzo+26 imrbi+10 iadz1i4+30 itsz2+25 
icspi Lavz2+284+e fadzi+28 
iet2or Aspo r3, ,1ta37r 
imrb1 ldvpi+e5 idvzi itst 
isch ier13 icat imrsi+2+c 
Ldvzo+ec ladet 1tsZ3+e+¢c 
Lmrt imrsi+#28 its2t icapi+2& 
4dvz2i4+20 imrg3+26 idvz24+26) iadot 
1eri0 {caz3+26+c Lav Z2+28+c 
iet33r ispO r+ ,ita75r 
imrpd1 fdavpi+25 i1dvzt itst 
fdvzie32 itset imrz4 fadxt 
isc4 icr13 icact imrzt+2+c¢ 





jadxi¢2¢+c imrt 
1tsZ3+2+¢ 
Amrt imrzi¢22 
1tst 4caz1+4+35 
f4adxi¢+28 imrt 
awe Ot icap1 
icabi410 imrpi+28 
jad4t iadot 
imrz44+30+c 
icaz3+28+4c 
1et65r ispO 
isco tera4 
tdvzo+46+e 
TELM | r6O,itaiir 
. cafl+rlt8+c 
coe aol 
et - E itsf34+48 
jadz34+46+¢ 
iede OT 5Y 
itsf4448 iscé 
sdvz2+46+c 
ispO ro,itai3r 
imrbi+ic itsxi+¢+c 
imrbi414 itsxi+30+c 
isce Lere25 
rie, ,itaeir 
icagitc  imrt 
jer10 icabi+10 
1dvet itszt+30+c 
fer15 
icasc+c itssit+c 
icaii+4 iadp1 
ri4,isc3 
ierid5 icase+c 
isupi idvpi 
icaii+e iadpi 
isupi+36 icpai 
idvpi+i8 imrpi+c 
ail,iscoO icr9 
{spre ispr3 
jeabi+18 imrfite 
icace+c itszitc 
itszZ14+3ce+c 
ier2 icafi-c 
ispr ispro 
imrfc+c fadsco+e 
icac3+c itszit+c 
itsz1+3e+e 
lere icafe+c 
ispr ispr7 
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Whirlwind Program 


its4t icaz3+c 
icti7r icact 
itshtt {cab1+10 
imrgi¢4O idvzi+34 
Its6t fceab1+10 
idvpi+2 iaddt 
lady idvzi+20 
itsz3+28 isc5 
Ladx1+30+c 
Ldvz2+28+c 
r5,itaiir 
icafi+4o+e 
itsfi+-46+c 
icapi+32 
fadz3+46+e 
gcaec LO6g5r 
ise icre4 
Ldvzo+46+ce 
ispO r3,italir 
icro4 icaft+48+e 
1tsfi+46+ce 
isc3 
1et3r isc5 
ict9r 
icapi+32 itsxi+c 
jicapi+28 idvbi 
ldvpite itszt+c 
imrpi+c jadbi 
1ecisr 
icter icait 
itsii+4  i1isup1+37 
itssi+ec icter 
jadpi imrb1 
itsiil+2 icaii+6 
ispsic #2,1isc3 
itssi+c 4eter 
icaci+c itszi+e 
ispr5 ispr9 
itsse+c ictyi+2 
ictae+e isct 
icta2+7  icaci+16 
itszt+c ictad+14 
isprio ye,isc3 
jtsse+e ictyeotre 
icta3+e ac 
icta3+7 fcaca+16 
itsz4+c icta3+14 
v3 ~44¢5 iert15 


Ldvzaet+e 
imrz4+28 
imrbi+210 
itset 
idvpi+10 
imr6t 
imrz3+26 
iert0 
imrt 
jaadt 
icapi+32e 


f§adz3+464c 


icde 
itsfo+48 


ispoO 


teaf3z+48+e 
itsf3+4G6+ce 


ied2 
icra5 
ierio 
ispo 
ict3r 
itst 
4ct6r 
itset 
ispO 


iadpt 
icpat 


icat 
itsbi 
o 


jad4tt 


iadx14+28 
idvzi¢e2 
imrzt428 
jar" 1 
its6t 
Ldavz2+26 
Loge: 
Ltsdt 
{ts73+30+c 
itsfi4+48 


hot5e 
ise6 


Lavzo+4 Gee 
P ialttal oe 


icapi+32 


jadz3+46+e 
5 Oar 
icafi+e 
ieafi+30+4+¢e 
r10,itacr 
ispO 

isc3 

is@5 
1caz4+30+c 
Pi 348 8 5 


icrt5 


itsit 
ispsil 


4MOA-123.1234¢ 
4MO04-123,1234¢ 
itsii+6 
icasi+c 


RCL 


wy pee 
a2c,isco 
icr 
its21+40 
ispri 
{eT 15 
a3, iscO0 
ior 
itsz14+40 
ispri 
icabi+4 


coe 
icr9 
ieacitdhec 
isce 
ispere 
icabi+e 
icr9 
icacl+4te 
ipec 
ispre 
imrf3+e 





jadse+ec 
itsZit-e 


itsZ1+432+c 


icre 
ispr+ 
ladsa+e 
s%,isc3 
icts9+e 
icapl 
jadet 
jadpi 
aOR 
icafitc 
itsfite 
ispric 
etic 
icazi+c 
2 bs Or 
isee 
{MOA-p1i23 
tothor 
ict4or 
iet5er 
Lits4t 
4etoor 
iMOA 
iMoA 

we LOL 

1 SOWRT 3 


+> 


itsse2+e 
icta4t+2 


icaf3+e 
—, 
itssc+ 
Pert 
#10, ees 
1adb14+16 
lepri3 
Leper 
[ 
itsgitee 
Youre 
ise. 
icaf3+c 
itsf3+e 
isprie 
Lone 
LOZ 4G 
isce 
16cc 
ical4 
isc? 
icait 
uaC3 
LOU ( oe 


ae 
ia 
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wnirlwind Program 


icty3+e 
isel 
ictat+7 
1tszt+ec 
yt, 1463 
Lety44o 
Leat 
Cm 
isuct 
ictsi04+2 
Tegel 

4-0 

eo Ger 
ise2 

i ogre 
itsg 
hemeGr 
1sc2 
icaltit+e 


iceafo+e 
icaf3+c 


lcaf4+c 
IMOA 


ores 


wt,iscO icrd 

rine icac3+4+c 
jeac3416 itezi440 
jetas+14 ispr2 
Ler15 Lcabi+6 
so,fi¢api idvs2418 
imps Qte itsse+e 
icasi+c Ldvse+c 
ema. ’3 icspi 
sil,icab1i¢+16 

Lsupi iadbi+16 
+25 sho is ec 
ispric i6Cce 
icaf2+e itszee 
icaz4+c itsf2+e 
Ict cur ispric 
isc? icart+c 
jeazgh+c itsfd+e 
idvf*. itsfie¢+c 
iavr4 ltteer ove 
jdvr4t Lusf3+c 
havAt itstt+e 
icalii+2 LMOA 
icas@tc = imrpi+18 


icack+e 


Lee 
ispre 
Jrrf +e 
1 Get 


ite set 
iadb1+10 
peut 
icpri4 
1 OmeES 
icaz4+ce 
cc i te 
eT 1 ae 
isce 
itszitec 
Lega oT 


1.MOA 
iMOA 


LMOA 
icapnl 
ilavpite 
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Table IV of Appendix B gives the Wnirlwind input data 
for runs (1) through (10). The Whirlwind program for run (8) 
has been given on the preceeding pages. The storage assignments 
for the input data are as follows (using run (8) as an 


example): 


Numerical 
FLAD Value Description of quantity 
pt, 60.0 R, assumed core radius (cm) 
0.052472 (v2,) 
20%.0 fits hefiector thickness (cm) 
20 Ar' = T/10 (cm) 
6.0 (X25. )y , inelastic scattering (core) 
0.0 a) , inelastic scattering (reflector) 
Cao3 . "{teration criteria for 3% accuracy 
0.000219 (»#,), 
f 
en, 22209 D, diffusion coefficient, core, Group I 
1059 D', aiff. coef., reflector 
0.621 f, f factor, core ’ 
0.68 f', f factor, reflector ‘ 
0.02158 fF) for core 
0.06c84 E45)” for reflector . 
0.000406 a9 for core absorption " 
0.0 =,’ absorption for reflector " 
e - =awal 17 Oa, reciprocal of Group I lethargy width 
ce, (See program) Same quantities as in cl, Group II 
03, (See Program) " . “= ve Group III 
ox, (See progran) - : — se Group T (thermal) 





() 


(2) 


(3) 


(4) 


(5) 


(6) 
(7) 


ae 
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CHAPTER II 
Numerical integration Method 
Re @. Kennédy 


I. Derivation of Equations 


The followins diffusion equation is assumed to 


zr 


apply: 
t 2 | 
3 ar, V Di = 2 abs; B - ~_ 082 5; Ai (1) 


3 1Qe ne 
The subscript i refers to a particular energy group. In 


this case (considering six energy groups) i ranges from 
1 to 6, where grouv 6 neutrons are of highest energy, 
group 1 of lowest (thermal). 

For groun 6 nquation (1) is modified to include 
a fission neutron source. All fission neutrons are 
assumed to be born in group 6. 

nauation (1) is to be made dimensionless, con- 
verted to a standard form (to allow the computer to use 
one standard calculation form for all groups) and integrated. 
foets first multaplied by (Lref/2ref) , where Z2ref is a 
reference macroscopic cros3-section, here taken as the 
thermal transport cross-section of the reflector. 


Bauation (1) becomes: 


“This derivation i3 based on a method developed by a. S. 
Thompson (J. Appl. Phys. 22, 1223) with modifications intro- 
quced to include leakace of the fast group (6) and in order 
to exvloit comnuter advantages. The orisinal method was 
developed for hand calculations. 





ee 


1 2 ref gy? Dabs 2ref ser. 2a, 
3 2ref Ztr, Yb. = Pref D5 - = reek = i gret gle! 


lore 


For convenience the follovine dimensionless cross-sections 


are definenr: 


a, = 2ref/2tr; 
b; = Zabs;/ Zref 
ds = 23; / mret 


2f5/ Sref 


ay) 
’ 


with these definitions haquation (2) becomes: 
a5 vy, ee ee _ OF di Zyros i (3) 
3 2ref 1 “rel “1 & log,b 
Next, to simplify the energy variable the follovine defi- 
nitions are e3tablished: 
6 e 10S, (Epi 33/k) 
1086 (heiss/Eth) 
¢ 6 (Los, Keiss/Een) 
eS, 8 
YAtogé = Ke x °Vd log E 
= 134 Vag gus 


eon © ag, 


anc the operator (-c 0/d@) may be 3ubstituted ror 


the onerator (dO/od log .f) in Mouation (3), 21V ime: 
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— 7 "By - 2 g - 2542 reths (L) 
aa a x rer he) 
Finally, the operator Thee which is given in spherical co- 
ordinates (with radial variations only) by 
me = o/s 2 ue & Sr 
becomes, upon considering a dimensionless space parameter 


r/R (where R is the total radius, including core, reflector, 


and extrapolation distance): 








2_ 0 Z 0 
Vv 3GAR * CAR OAR 
(ey 2 2 d 








ae 
"OC  “E@7R) 8672) 
Thus Equation (4) becomes 


a 2,2 O& ca,2.,. 
(R°V “)f, - b, 2.8, 2 CS “ener 
Saar 4 1 “ref 4 ee (5) 


or, Since the operator (R°Vv*) In rapace is equivalent to 
the * ee “a in r/Respace: 
3 a RS ec: ie 
Upon multiplying Equation (6) through by 3( 2 BR)» and 
defining 

i Ee oh Pret)” 


Equation (6) becomes: 


cd ( 
a,V* f,-»b Ae, = Ob AotaAa 7) 
Note that the parameter Av depends on the total radius of 


the reactor and is thus not originally known. Instead it is 


assumed (based on one-group calculations) and later corrected 
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ii the assumed radiu3 is subd- or super-critical. 
fhe derivative on the right hand side of Equation 


(7) may dé aprroxinated: 
DRE CAH, ABCE EO “ABCs Ke 
@ A 
Boone aiVi “bs a by AY; a Bcd, & = neCde, ea (33) 
AG Ae 


The left hand side of pguation (8) represents loss fron the 
1*th enerty sroup by Leakage, absorption, aril slowing” down, 
Poe rient and Side represents the input to the i'th grow 
consisting of neutrons Slowed from the next higher croup. 

Nor i=l (thermal group) Equation (@) 15 modified 
to: a, Vp, ~ b, A 2s = Ae cd 5 (8a) 
Since no neutrons are lost from the thermal grouo bv slowing? 
down. 

The input term in hquation (8) is modified for 


eneresy groun 6 (fast) to include fission: ; 

ap 7b, - be AB = we B cand, = AVerch) , Davi iYs (9) 

ae ae La 

(The fission of thermal neutrons must be weighted by an 
energy factor since all other fluxes in (9) 2re per unit 
enerev interval and must be multiplied by this encrazy inter- 
Pet to become tha total flux. rhe equation is Sarolittegar. 
dividine through by the enersy interval with the result that 
it apnears as a denominator in the therpel term.) 


pquations (&), (€a), and (°) are of the rorm 


ie rs R = 
(a.V = te Bs =o h 
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where ss combines all lo3s terms except leakaze, and h repre- 
sents an in»out term. ror example, considering energy zroup 1, 
e@ is siven by (b, A) and h by (=A zo) from (8a). 
Ir this form of hquations (8), (8a), and (9) is 
Multiplied tnrough by a.(r/ay” and the operator av" 
replaced by its eauivalent (a; 1 d/d(r/n) (Ji)? A/a 
the tollowing expression (vith cor Bhi results: 
a,x" {(1/x2) a/ax(ayx? )d/ax - g* } J, = a xth 
or, with v, = asx" d7idex, 
i: - a; he ) ) Js = ext (1U) 
Lquation (1U) is of the tvpe 
(De - axe”) Y = 2 


which may be factored into 


(poe f(D = £) YSZ (11) 
provided f satisfies 
af ies 
Df +f - axe = 0 (12) 
This arisés from 
(D4 f(D -f) Y2zs2 (De = axte%) Y 


p2y + epy - D(fy) - £%y = Dey - ax,7y 
D*y + py ~ (DF)Y - £(nyY) - £2¥ = D°Y - ax 
or, (DE)Y + (£°)¥ = (axde~)y = 0 
me Zives houution (12) specifving f. 

In heuation (11) let (D - f)Y = U, so that 


(pear )U = 2 (13) 
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where (DD. Be. $8 (14) 
Integration of Equation (14) gives 
xX x “a ‘ ' 
Yo exp ( 0 a. axp axe’) { f exp(- /*F aR aE) U a + e } 
(15) 


and from Equation (13) 


U= expl- ABR SRY (enol (Eats) 2) aatts 5 2 c, } (16) 
© ax! aX Ea aati” 
which combine to give (writing m= - f ax" /ax'*) 


Yo=n fFor*) ( m Z dx! ¥/axt#® + Co] axt! t/axtt te + Cy 
(17) 


The constants C,; and Ca are determined by the conditions 


that 
y' x20 atxe-= 0 and 


Y= 6 A@ xe 1 
which is equivalent to requiring neutron flux symmetry at 
the reactor center, and requiring the flux to vanish at 
the extrapolated radius. 

By the first condition Co must be zero. By the 


second ©, is given by 


5 * i (1/m*) Vo Z ax! /axt® Jax /axt? 


so that Equation (17) becomes 


Y= of f*|(a/n? a: Z axtfaxt® [axt¥/axt'?} (18) 


Solution of Equation (18) requires prior calcue 
lation of f from Equation (12) with the initial condition 
that f(0) = 0 from (11) and (10). Equation (12) 1s written 








in terms of differences: 
‘ 
ax*d/ax(f) * (f)* - axto*~ = 0 
or, with subscript 1 indicating space point x and subscript 


Patndi cating (en2ce oin’ (x—-— Ax): 


es 
ea nes os Za 


2 
(fo fo x = a(mo KD + GIy ) - 2(2>5 + fae 
Solving by quadratic formula gives 
f, = Hanx/b x + apt, { a 5305” fen) © + pty, Oe + oN Bg + 
Ze shies te $ 
+ FX" = Py" /ayxy \ (19) 


Factors in squations (18) and (10) are summarized tor 


elaricy: 
Energy Group be (ameut ) 2° (loss) 
1 (th) (~A 3 ed.) b A 
2 (-A® cd, ), boA 4+ (cd./ae JA 
AR 
3 (-A& cay ) 3), b,A rs (cd3/ ae AS 
4G 
L (- AS cd. )B, br + (cd, /ag rX 
f 6 
6 < ¥ e10Py - ve beA + (cde/aP rE 
AG 


1 
Il. Froerammineg of Problem for Computer 


Whirlwind I (WWI) is a digital computer canable of 


performing simple mathematical overations (addition, sub- 
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traction, multiplication, division) with extreme rapidity 
and under the control of a set of instructions, called the 
program, which may be stored in the computer. WWI normally 
operates with numbers that are all less than unity, but, by 
accepting a loss in computer speed and storage capacity, may 
be operated with unscaled nunbers (thereby decreasing 
enormously the difficulty of setting up the program) in the 
Comprehensive System (CS system), as ig done in this problem. 

In the CS mode there are 1360 registers available 
for storage of either instructions or numberse an ine 
struction occupies one register, but a number occupies twoe 
Of the two registers containing a number, the first carries 
a part of the significant figures which represent the 
number expressed in exponential form, while the second car- 
ries the exponent and remaining significant figures. For 
example, the number 123.4 may be expressed as .1234 x 10°, 
with .123 stored in the first register, and the exponent 3 
and the remaining significant .0004 stored in the second. 
Actually, WWI operates with binary numbers, but an under- 
standing of this feature is not essential to the under~ 
standing of the problem since the computer accepts decimal 
numbers from the programmer and returns decimal number 
answerse 

The 1360 registers which comprise the internal 
storage capacity of WWI in the CS mode are called the 


magnetic core memory (MCM). External storage is available 








[ee 


in the form of magnetic drums and tape but acquisition time 


is large in comparison with the time required for the com- 


puter to obtain ae number stored in FCM and is therefore 


not used in this problen. 


Any problem then is limited to a set of numbers 


and instructions which occupv not more than 136U refisters. 


In this problem these registers are used (starting at the 


first register available, which is number 32) as follows: 


Regzisters 


registers 


Registers 


Registers 


registers 


Kheristers 


32-209; 


21-460: 


L70—-5093 


510~597: 
5086630: 


640-end: 


working storage (for numbers 
which are calculated during 
the problem and must be ten- 
porarilv stored) 

storage for 20 values of the 
neutron flux in each of the 
six energy groups 

Storage for 2U space points 
(r/x @ -05, .1lu, .15 etc.) 
cross-section storage 
Storage of constants other 
than cross-sections 

Storage for the program 


instructions 


Actually the program for this rroblem uses onlv some 1¢5U 


registers so that 110 are available for later expansion of 


the problem if desir 


ed, 


For the purposes of this problem there are 31xteen 
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commands which wWwI obeys in the CS mode. They each affect 
an operation, designated by three letters, and an address. 
The address may consist simply of the number of a register 
which contains, for example, the number which the computer 
is to obtain from MCM.e However, since a programmer cannot 
know in advance exactly how many reristers to set aside, 
it is frequently difficult to specify a register number as 
an address. To altace this problem, WwWl 1s designed to 
allow use of indefinite addresses, called flads (for 
floating addresses) and represented by a letter-number com 
bination such as al, a2, bl, b2, ete. Thus, if a number is 
to be stored at the end of the program it is tagred with a 
flad and the flad used in place of the actual register 
number. WWI simply counts the registers as the program is 
read in and remembers the actual address of the flad. fThen, 
whenever this flad is referred to later Wwl substitutes the 
proper addresSe 

The three-letter portion of the WWI commands gener— 
ally calls for a certain arithmetic operation. All of these 
oceur in the Multiple Register Accumulator (MRA). For ex- 
aap, if a number in storage is to be multiplied by another 
stored nutiber, the first is taken from MOM to the MRAe The 
second is then obtained, the required product formed in the 
MRA and then transferred to storazee Various operations, as 
explained below, may either preserve or destroy the originally 


stored quantities. 





ah 


The sixteen instructions used in this problem are 


the following: 


Instruction Meaning O ot duon 
its al transfer the cone the number in the MRA is 
tents of MRA to transferred to register 
register al al; the former contents of 
al are lost; the contents 
of MAA are unchanged 
lex al exchange the cone the contents of the MRA 
tents of MRA and become the new contents of 
register al register al, and vice versa 
fea al clear MRA and add the MRA is cleared (and 
to it the contents its contents lost) and its 
of register al new content is the nunber 
stored in al; al remains 
unchanged 
ies al Clear MRA and gub= as above except that MRA 
tract from it the acquires the negative of 
contents of al al; al is unchanged 
iad al add the contents the contents of al, plus 
of al to MRA the original contents of 
MRA are added and the 
result stored in MRA; al 
is unchanged 
isu al subtract the con- the contents of al are 
tents of al from subtracted from MRA and 
MRA the result stored in MRA; 
al is unchanged | 
imr al multiply MRA by al the number in MRA is mule 
and round off tiplied by the number in 
al; the result rounded 
off to the maximum number 
of significant figures 
which can be stored (about 
8 decimally) is retained 
in MRAs al is unchanged 
idv al Givide MRA by MRA is divided by con- 


contents of al 


tents of al, with the 
result stored in MRA; al 
1s unchanged 
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The following two instructions provide transfer of 
control to the instruction stored in the address of the 
instruction. WWI normally operates sequentially from one 
resister to the next. These two instructions circumvent 
this sequential operation. As an example, the last ine 
struction in an iterative problem might contain a transfer 


of control back to the first in order to perform the 





iteratione 

Instruction Meaning Oneration 

isp al control transfer take the next instruction 
from register al and pro~ 
ceed from there 

Lep ad conditional as above, if the contenis 


control transfer of MRA are negatives other- 
wise, ignore, and take the 
next instruction in sequence 
The next five instructions have to do with operations 

performed with a counter. In WWI a counter may be set aside 

to allow repetition of an instruction, or set of instructions, 

a specified number of times. The svecified number is the 

criterion, the number of times the operation has been per= 

formed is the index. WW1l sets aside both criterion and 

index (with index initially zero), cycles the index once 

each time the operation is repeated until, upon comparison 

the index is found to equal the criterion (or if the index 

exceeds the criterion, as pointed out below) at which time the 


computer proceeds to the next instruction in sequence. 





lige 


nearmuction Meaning Overation 
ice j select counter j select counter j for fol- 


lowing operations with 
Counter 


(selecting a counter J, by isc 3, does not in itself set 
criterion and index; this must be done by an addittonal 


step prior to the first vee of counter j.) 


instruction Meaning Operation 
her mi cye le» reset set criterion equal to m, 


index zero 

ict al cycle count increase the index register 
by one and compare with 
Griterion, ltetme index is 
begs return to register al: 
otherwise take the next ine 
struction in sequence 

fca fg cycle increase incmeasie the index by a 
(this may set index greater 
than criterion, but is use= 
ful in operating with more 
than one counter) 

een cyele decrease decrease index by n 

(Note that the counting instruction, ict al, is the only 

counter ecmmand that transfers control. It is therefore 

usually the last instruction in an overation which isa to be 

repeatad.) 

During operation with a counter WWI may treat 
numbers (or, in special cases, instructions) which are stored 
in sequence if +c 15 appended to the address of the first 
number to be treated. For example, to place ten numbers 


which are stored consecutively in ten other consecutive 


registers it suffices to select a counter with criterion ten, 
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@adex zero, obtain and transfer the first number, with ¢c 
apoended to the address section of the obtain and transfer 
Paotructions,sand transfer control back to the first instruc- 
tion (in, say, register by) with an ict bl. hach time the 
coInouter notes the +c in the address it sets c equal to 
twice the index of the counter in use and thus transfers 
the ten desired numbers. The commuter sets c to twice the 
index hecause of the fact that the numbers occupy two 
registers. This however is not the case if a +c is added 
t) the address section of an isp command. Here the address 
mMuSt of another inetruction, occupyine one register, so 
the computer conveniently sets c equal to the initex. 

The final instruction is used in connection with 
subroutines within the orogram. Identical operations which 
are to be performed at different points in tne propram nay 
be Set ud aS a subroutine in order to avoid duplicition of 
instruction storage each tine the recuired c.ilculation 1s 
Wem. hen, at each point the sudroutine is needed, an isp 
instruction transfers control to the subroutine. Howcver, 
to leave the subroutine and return to the proper point in 
the muin program an isp instruction at the end of the subd- 
routine must have the address of the rezister inmediately 
folleywane the r@gister which originally transrt@rred control 
to the 3ubroutine. whis is accomplished bv the rinal instruc- 


PpiUemeEO be cOnclacre. here. 
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Instruction Meaning Oneration 
ita al Sransfer address the address section of the 


instruction in al is changed 
to the register immediately 
after the register con 
taining the latest isp ine 
struction which has been 
obeyed 
For example, to obtain the square root of a number at several 
points in the main program, a square root subroutine may be 
set up at the end of the program. Then, each time the 
operation is to be performed the number is placed in the 
MRA, control transferred to the subroutine, the root calcu 
lated in the subroutine and placed in the MRA, and control 
transferred back to the main program to the register ime 
mediately following the original transfer point. 

In the subroutine the first instruction would be 
ita bl, with bl the last register in the subroutine. In bl 
would be stored isp -—-. Thus the ita instruction would 
cause transfer of control back to the main program at the 
proper pointe 
ITI. Prorramning of Specific Problem. 

The problem to be solved here, derived in Part I 
of this chapter, may be stated as 


Y¥(x) = rae | = ) exp(—2 “f dxh AM | ox 
O at awe 


ax! 0 wu” 





Zax } Iv 
ax 


tw 


2 
ax tit 


where x 19 a dimensionless space parameter r/R total 








ee 


@is a 8tored dimensionless cros3-section; 


4 is the source term, representing fission neutrons 


(for the highest enerey group) or slowed neutrons 


(for all lower energy groups); and 


f is calculated from e« finite difference approximation. 


Fundamentally the computer proceeds by the fol- 


loving &tens: 


ag 
(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(&) 


Calculate f(x); 
Nivide f(x) by ax* and interrate lrom U to 
x by trarnezgoidal summation; 
- 2 

Calculate and store exp ia ae 

O 
Calculate Source term 2, initially from 
an aSsumed flux and critical radius; 


£ and multiply by the expo- 


Divide “4 by ax 
nential computed in step (3), preserving the 
stored exponential for later rr -use; 
x 2 Xx 2 
Integrate (if/ax~) (exp t dx/ax ) dis 
O 
NYivide the above interfral by the suuare or 
the eryponential (still preserving the expo- 
Bec P ne 
nentiel) and by ax”; 
neverse the order of storace of tne Guantity 
determined in sten (7) so that the intescral 
fron 1 to x mav be determined (this step 


13 necessary since t1.I counters increase 
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rather than decrease following each repetition). 

Step (8) effectively cancels the minus sign which 
appears in the input term to each group, with the result 
that the calculated neutron flux 18s positive, as it should be. 

(9) Reverse the order of storage of the integrals 
determined in step (8) and multiply by the exponential term. 

(10) Store the results of step (9) and proceed to 
the next energy froup calculations, changing f, a, and Zin 
steps (1) through (9). 

Following the sixth application of these steps the 
new thermal flux has been computed; wWwWI proceeds then to the 
next step: 


(11) Compute a new critical radius parameter A!: 


’ ‘ J 
To , g,,(xiar/ [ 9, (ne A 


where the primed quantities are the newly computed valueSse 
Step (11) is based on the fact that once the 
proper solution for the flux is determined D., will not 
change with more iterations, as long as the radius is equal 
to Roerpiste On the other hand, if the radius assumed was sube 


or supercritical Ben will be less than or more than Ze, 
(assumed) respectively. 

(12) Compare the two values of the radius parameter; 
if they differ by more than an acceptable error recalculate, 
based on the new values of flux and radius; if the reactor 
is critical, as evidenced by a specified small variation in ma 


-t 








wm Si 
print out specified information (e.2., flux, radius, power 
distribution) and stop. 

The computer mav be stonped at the end oft a 
problem by an iSTOP instruction. Alternatively, if more 
problems are to be considered they may be listed after 
the first problem and the comvouter stonped with the special 
instruction so26. when it reaches this command at the end 
of the first problem the comouter automatically reads in 
the remaining problems consecutively. in this way a 
considerable saving in time may be realized, since manual 
effort (otherwise required to start each problem} is mini- 
mized and since many of the constants which are already in 
Storage will not change from problem to problem, thus 
decreasing the time required for the computer to convert 
the innut decimal numbers into binary numbers. After the 
final nroblem an i35TUP command may be read in to stono the 
comruter. Utherwise the overator will stop the computer 
manually when, upon ccmpletion of the last problem, the 
computer, under the sp26 instruction, calls for more infor- 
Pies 2 ON). 

A3 Mentioned in part JI of this charter, the 
numerical values of this problem are stored in the first 
toon WL OMe. The first inStruction is stored in register 
GLO. Develupment of the srogram for this preblem will start 


at this point with numbers specified as they appear. 





Instruction 


azi,iscO 
tcrept 
isel 
lerppe 
isce 
Lered 
ise3 
ier2o 
isc4 
icro 


1¢ca636 
its598 


igc2 
a3,ica430+c 
imr63 
its430+c 
1CCA 3 


f2,isc0O 
a2,icad510 
imr254+e 
its34+c 
ica51e 
imr298+c 
iad34+e 
its34+c 
icad514 
imr342+c 
iad34+c 
its34+e 
iea515 
imr386+c 
jad34+ce 
its34+e 
ica518 
imr430+¢ 
jad34+e 
imr598 
imr470+¢ 
imr470+c 
its 34+e 
icvac 


5 5 
Operation 


Sets counters: O for iteration over 
core region, | over reflector, 2 and 

3 over entire reactor (20 space points) 
and 4 over six eneryy groups: ppl and 
pp< are preset parameters which are 
set at the op ning of the program in 
order that the rest or the program 

may remain unchanged 


Stores assumed radius parameter in 598; 
636 contains latest calculated A 


3cales Ben to any desired center-line 
value (value initially set in 
register 638) 


Selects counter U to iterite over core; 
double-length registers 510 through 518 
contain fission parameters; 25/4 through 
L6& contain fluxes; input term Z to 
fast energy group is determined as the 
sum of each (y B; ) miltiplied by 

(A) in 598, x©"(x in 470 through 508) 








ica574 
its618 
ica576 
its620 
ica530 
its622 
ica53e 
its624 


ispmt 


a+, isc 
ictm1L 


ispx6 


md ,iacO 
a5,ica618 
imr470+¢ 
imr470+c 
its78+e 
duc Ca 


isci 
a6,ica620 
imr470+pp3+c 
imr470+pp3+ce 
its /Crppjt+c 
icta6o 


ica620 
iad618 
imr614 
imr470+pp3 
imr470+pp3 
its78+pp3 


iscO 
a7,ica6e2 
imr470+c 
imr470+c 
imr59 
itsieet+e 
1e¢ba'7 


eed 
a8,ica6be4 
imr47O0+pp3+¢ 
imr470+pp3+¢e 
imr598 
itsil22+pp3+c 
icta8 
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Stores loss parameters (b+B cd/4®) and 
a (2rei/2tr) or group 6 in working reci- 
sters 618-621, (this is done in order 
that other group paruneters may later 

be 1.laced in these sunie working resi- 
gters and the remainder of the program 
used for all energy gfrouns) 


causes a jump around the next 3 instruc- 
tions which provide entry for groups 1-5 
ay is ihad to lJadel entry point rom crops 
1-5; next register jumped unless new a 
has been culculated 

jumps to %,, subroutine in x6 when new 
Yi. has been comvuted 


2 


ealculates ax~ for corq region 


as* above, for reflector refion; pp3 is 
set at the opening or the program to 
twice (because of double-length reyi- 
ster3s, the estimated value of Rtcore/R, ., 3 
thus, in order towpecizy « refleetor 
core combination it is only necessbury 

BO changes ppl ,<,3, @9ol progral 


a 


places av rug: vulue of ax~ at core- 


reflector interface: 61). contains 0.5, 
used frequently in averaging 


caleulates. loss tern (.e%., bA + Acc /a8 ) 
times x* for core 


ag above, for reilector 








1¢5122+pp3 
ica626 
Ltas 
LEB EOC 








eed 
jer’ 9 | 
ab acsi6Grtc 
imr16o+e 
imr470 

Lav 734e 


1adi6&+e 
iad16¢+e 
imr+7C 
i. dbo+-e 
idv47o 
idv4790 
Jadl2ise 
jad1l2e2+c 
4mr80+e 
ispze 
itsi70+c 
LesS0+c 
idv470 
1adi7C+e 
Ltsi7%+c 
ictag 


@.,ical6l+c 
idv7G+c 
itsi6c04c 

1 ome 1 


ae Lace 3 £4) O in vegisters 325 129 
oe in miLly se— to zero and ler 


to @et as @ero gource 


-~ 


SBloulatvey fT at xe. 35> (Met be calculated 
SUpemase Ly to avoid jivisio: by yen 


ae 1c. SQuUaPS reOt. MUupresa line 


calculates remaining 19 values of f 





calculates t/ax® 








isc2 
b1,icaiLSo+c 
{adi68+e 
imr6o41 4 
imr470 
itslee+c 
Level 

isce 
beyicale2t+c 
iad120+c 
itsil2etc 
ictbe 


isce 

be, cad @2+e 
isp2i 

i medcct+c 
Lee 3 


1Ssec 
b5,ica34+c 
imrieet+c 
itsa34+c 
Lemb5 


isec 
bO,ica32tec 
iad34+c 
imro44t 
imr470 
its168+c 
ictb6 


ica626 
1ts166 


LSee 
b7,icaL68+c 
iad166+c 
its168+c 
meu 


iscO 
fi,ical68+e 
idvieet+c 
idvLleet+c 
Ldv#70-+c 
idv470+c 
Lav618& 
itsi68+c 
ietrt 


Je. 


forms areas of f dx/ax* Prior te 
trapezoidal intageition (1.70 =A x 


sums above areas in trapezoidal inte- 
gration (120 previouslv set to zero) 


x 
2 : 
calculates oot [ fide ei lag 
21 identifies O 

exnonential subroutine 


source term Z times exponential 


forms areas of (Z times exponcntial) 
prior to trapezoidal inteszration 


3ets 166 to zero 


Suis areas, in trapezoidal integration, 
to compute fx 5 
Z (exp) dx/ax 
O 


for core region: 


Z (exn) dx/ax* 
O 


vs a 
(ex) .ax 
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isct 
@3,1ca168+pp 3+c 
Ldvi2e+pp3+c 
Ldvi2e+pp3+c 
Ldav470+pp3+c 
idv470+pp3+c 
idv620 
its168+pp3+c 
aete3 


{ca620 
4ad618 
imr61 4 
iex74 
Ldvi2e+pp3 
1dv1l2e2+pp 3 
idv4+70+pp3 
idv470+pp3 
idv7it 
1ts168+pp3 


isc3 
icr20 
icii9 
ke ise3 
icai68+c 
Ledd 
16Cc2 
its34+e 
ictkt 


4ca626 
its 7/4 


isee 
b8,ica36+c 
iad34+e 
imro14 
imr470 
its/7ore 
ictbs 


1ca626 
ivs'*6 


2585 


Je 
sets aside (in my f i (exp) dx /ax” 
> 


at core-reflector interface 
for reflector region: 


é 2 
Z (exp) dx/ax 
b 


(exp) * 


aS above, at core-rceflector interface 
COmLNe numerator Sate asideein 7/1} 


reverses storage of above (over entire 
PesctOsy Tigo de™ bo interrate tie 
libs vig mee: 


Sets rnecister // to zero 


forms areas prior to integration 


sets register 76 to zero 
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isce sums above areas in trapezoidal inte- 
b9,ica78+tec gration 

fad7b+c 

its78+c 

ictb9 





Ese3 returns above intesrals to original 
icr20 Omeer ot storage 
icit9 

k2,isc3 

ica7O+c 

icdd. 

isec 

its34+e 

we tHe 

isce 

k3,ica34+c 

its78+c 

ictk3 


isce2 commutes flux: 


ct ,ica7St+c x ia Ligaen 
imr122+¢c GZ = (exp) Z (exp) dx/ax E 
° ey 

1 i ,  bexp)® ax” 


isc4 counter 4 selects energy group by jump 
isppl+c Go xl three o> ivnech store Be through 
pi,ispxt J> and position constants for next com- 
Lspx2 nuting Ye through J); the isp ab 

ispx3 instruction resets counter 4 and pro- 
ispx4 ceeds to the J, storage routine 


Lsepxr5 

ispa4 

x6,isc2 to determine radius criticality calcu- 
10a626 lates an apnroximaution of i B(x) dy. 


c2,iad34+c 
tet ee based on the nev ly connuted Gy in the 


itsii6 recistem® following 31 


b. 


isc2 2S above, based on praviouslv computed 
‘aaa (or assumed) Jy 


c4,iad430+c 


ietc4 . 
iav1i16 computes ratio of above intcrrals ¥ 


isu508 comoutes YF- 1 

imr--- multivlies (¥- 1) by any desired tactor 
tad508 to Control rapidity of Cconversien ,~pius 
imr598 COMPUTES New 


1ts636 A, =, i+f- Yn) 











ica34 
Lsuje. 


iad34 
1tsice 


iaee 

co ,ica34+c 
idvi22 
imrd3e 
Ltst32G+e 
PETC 


12a 30 

1dv612 

tepz2 

Ldvdoo0 
4MOA+123.456¢ 


idv536 
{adoo4 
lepat 


LOR 
+5 
+0 
+20 


BSCe 

Cf ,lcat30+e0 
LMOAF1 2345076 
ove | 


BSt 2 
03,1030 +c 
MOA+1 .2345f 
4 cite 


= 


linear extra: olated value of of 

: ag 
ay remevor @@iiter to normalize 
Corl Cire Yen 


nurmaLiae calcstated Y | to corres- 
pond tooriginaily chusent value (in 
pepaster 30) at center 


store normalized ¥,, in assigned 
reristers 


calculate Heotal from A 


print out Ryo4qz (Sample number given 
to specify number of significant 
figures desired} 

determine criticality bv comparing 
the per centage change in 

with an acceptable error stored in 
rezister 604; if error is too large-- 
elther plus or minus--return to al 
(the original starting point of the 
pro¢cram) and recalculate; if error is 
accentable, continue to next sequence 
of Inateae tons 


svecifies format of solution: blocks 
of 2U numbers arranged 5 per line, 
with a tah between numbers 


prints out 2U values of the calculated 
flux for each of 6 energy sroups con- 
Sidered 
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etCc10 


ise 
cit, lea2dsi+c 
Lear 1. 


isce 
ei2?,ica2ice+c 
LMOH41.2 345° 
HeoC iL? 


isco 

mi, 104520 
mrt jo+c 
Lte3z4+e 
ioad2e 
Imr38ot¢c 
4 Gas ec 
Loo tre 
iead2't 
Lmr3 ‘hdc 
fad3zlt+ee 
its34+c 
icad52eo 
imr2jorc 
fad 34+c 
ite3t+c 
ica528 
imr25't+e 
jad34e+c 
its 34+c 
ic tmfs 


ica34 
18u36 
{ad 3'} 
itstec 


calculate the conribution to the 

total power distribution from each 
proup (fission cross-sections-- 
weithted by A® /c for ., Ac ie. 57 









stored in 520, 538ay 526 528) 


linear extrapolation to determine 
DOWEr av Celwer 
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my ,icaz++ec 
idvie2e2 
IMOA+1 23 'F9OF 
Loom? 


=_ -) eg a 


4,18 52 
d1,icastte 
ese1Orec 
eco 
ic@s/S 
ivee13 
105530 
1t5620 
4009534 
1ts622 
10a536 
1ts62'} 
isco 
db ,ica3zt+e 
imrl7O+c 
imr'+7O+c 
Lmr598 
Lore 54 
Lts3te+c 
ictdo 
{ca3z++pp5 
itecO0e 
ec 1 
e2,icutsterppjtre 
tine 4 7O0+pp jre 
imr'}/O-+po ote 
imr593 
Lmr556 
LES3"srpp se 
ievéec 
103554 
paadd55e 
imroL 
Lmrsoe2 
Lari 7/0+po3 
imr470+pp 3 
ime 0 
its3+pp? 
Lspac 
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print out power distripoution 
normalized to unity at center 


insert stopping command (e.f., 
isTOP or sp26, as discussed 
On page 5h 


subroutine to be used after 
calculation of Bo: stores 


in assigned registers, places the 
parameters of the next lower 
energy group in the working 
reristers and calculates the 
source term Z for the next 

lover group; reenters prosram 

at entry point set aside on 

page 56 


(THE FOLLOWING PAGE CONTAINS 
FUUR SIMILAR SUBROUTINES “HICH 
PERFOnE THE SAKE FUNCTION AFTER 
UnLCULATION OF $e, $i, ¥5%.) 
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KGa, Mec O XBgeieue 2 wt pdlisic 2 x5,1isce 
d2,icaztt+. d3, Leagit+c a4, ica34+c d5,ich 34+c 
Lts254+c its2e)3+c its342+c Lts350+e 
ictde2 Teed 3 ictd4 Teo 
Lew582 ica538c {¢a530 ica594 
Lts61é Lts616 smersion Ra its618 
ica5o4 ica538 icasge ica5jo 
Ltso20 its620 168620 iteo20 
10a53a ica5he ica546 1¢a550 
itso22 its622 beb@c Wie 22 
ica540o Lead544 1ca5%8 $e2552 
itso24 itso24 itso24 C8024 

iscO LSCO isco i1scO 

a7, lLea3z4+e dB,icar+c I}, lcaz4+c ed. yLeas++c 
Lrr4t /O+¢ Lar+70+c imr470+c¢ imr:-7O+e 
Lart70+e far470+c imr470+c imr470+c 
imr593 imr598 imr590 imr59d 
imr558 imr562 imr566 imr5FO 
Lts3%#+ec its 34+ec {ts34+e its34+c 
iced ic ta ictd9 ictet 
ica34+pp3 ica34+pp3 ica34+pp3 ica34+pp3 
its602 18602 1ts602 its602 

lect isel ise ise 
e3,ica3t++pp3+c et,ica34+pp34+ce e5,ica34+pp3¢c eb ,icaist+pp3+e 
imr470+0p3+c imr470+o0p3+¢ imr4+70+pp3+¢ imr4+7C0+pp 3+e 
imr470+003+c Lmr4794+ 90 3+¢ imr470+pp 3+c imr470+pp 3+c 
imr593 imr595 imr598 imr593 
imr560 imr564 imr568 Ln ee 
its34+pp3+e Lts344+pp3+c Lts34+op3+c its34+pp3+c 
icte3 icte4 ictes5 icted 
ica55o ica562 ica5o06 ica570 
Lad550 iad564 iad568 1ad572 
imro14 ere imr614 imro14 
imr602 imr602 imroo02 imro60e 
imr'}70+pp 3 imr470+pp3 imr470+pp3 imr470+pp 3 
Lmr470+p 53 imr470+pp3 imr470+po3 imr4'/U+pp 3 
imr593 imr598 tmr598 imr598 
its34+4pp 3 its34++pp3 Lts34+epv3 Lts 34+ p3 
ispa4 ispa4 isoa4 ispa4 


4, exponential subroutine (LSR FU1) 

22, square root subroutine (LSR FU 2b) 

~ START av 640 (terminates read-in and instructs computer 
to take first instruction from register 640) 
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Sore 
FLOW SHEET FOR SIX GROUPS 


(Superscript i denotes energy group, from 1 to 6) 


Subscript n denotes space point, from 1 to 20) 






Posu tion Group 1 
Constants 









S Cdax ps3 


ee W% Q(x) dx 
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FLOw SEEST *OR SIX GROLFS (Continued) 









wormalize 





Paleo 









print 
Power 
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FLOW SHEET FOR SIX GROUPS (Continued) 


SUBROUTINES 








ae 


Cin tan Ii 


ASSULTS 


I. Regetor Comnositions Studied 


ae Core Materialse In order to simplify the 





problem, a ehcice was made 2s to the type of fuel element, 
coolant and moderator. The choice was based on a considered 
adaptability to high output mobile reactors. -.High power re- 
quirements, counled with a need for light weight and small 
size, dictated the use of hiehly enriched uranium and clad— 
ding materials which would withstand high, temperatures. 

The fuel element was a dispersal element of 
uranium dioxide in a matrix of stainless steel. Some calcu 
lations have been made on this tyne of fuel element, and food 
high temrerature properties are envisioned. The amount of 
uranium dioxide in the steel matrix is kept low in order to 
minimize radiation damage and to improve heat transfer 
properties. <A similar fuel element, with the steel replaced 
by aluminum, was operated at the Geneva conference on the 
peaceful uses of atomic energy. 

The uranium considered is enriched to 90% Um235. 
The density” of UO, is 10.95 gm/cce 

Stainless steel is also used - the cladiing for 
the fuel element. Its choice is dictated by the comer that 


it will withstand sodium up to about 600 F. Ite absorption 





cross section presents little problem with the fuel enrich- 
ment considered here. 

Berylliwa is chosen as moderator because of its 
good moderating properties, light weight, and cood thermal 
properties. It is protected from the sodium coolant by stain- 
less steel, since beryllium witnstands sodium only to binauie 
1000°r and higher temperatures tnan this are considered. 

Choice of sodiui as the coolant is based on its good 
neat removal characteristics and high boiling point. A nigh 
boiling point eliminates the need to pressurize the core, 

Tne reflector is considered to be cozposed of 
various proportions of Be, Na, and stainless steel. Reflector 
cooling is necessary since some five per Bete of tne heat 
may be liberated in the reflector. Tnis amount of neat will 
De considerable for hign outpwe reactors. Steinless steel 
is required, as in the core, to protect tne oeryliium from 
tne corrosive ection of high temperature sodium, and for 
structural oDurposes. 

With a basic reactor cnosen, variations in the 
proportions of tne materials were studied in order to 
determine their ei fect on tne power distribution. Table III 
of ADDendix BE summarizes reactor compositions studied. 

The effect of reactor poisons was not considered. 
nl. of the reactors were assumed to be ho ogeneous and in 


the steady state. 
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b. Cross Sections of Reactor Materials. The 





cross sections of the materials considered are listed in 
Tables I and II of Appendix B. They are based on data 
from “Neutron Cross Sections," BNL 325 (August, 1955, 
Geneva Supplement). 

c. Reactor Configuration. For reasons of mathe~ | 
matical simplicity the reactor considered is spherical, 
with a reflector of uniform thickness. Because of angular 
symmetry the spherical reactor may be resolved to a one- 
dimensional problem. Cylindrical or rectangular geometry 
would, on the other hand, require twoe or three-dimensional 
analvsise While it may be argued that spherical reactors 
will not find extensive use, the results from an analysis of 
spherical reactors can be extrapolated to cover reactors of 
other configurations. 

The magnitude of the multi=dimensional problem 
is such that simplification to a single dimensicn is made 
frequently. Other than the spherical case considered here, 
tne only one-dimensional geometry is that of the infinite 
slab. Spherical geometry is superior, since it gives 
results which may be applied to reactors of other configu 
rations whose dimensions along various axes are not greatly 
different. 

de Fuel Concentratione The proportion of fuel in 
the reactor has an important effect on the radius of the 


critical reactor. It is also related closely to the power 
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level. 

For the mobile power reactor envisioned here, 
there are certain restrictions on the size which must be 
considered. t cannot be go large that its mobility is 
questionable, and yet it must be large enough so that 
enough coolant may be pumped through to remove the required 
heat. 

From elementary heat transfer considerations 
(Appendix A), a core radius of from two to four feet was 
found to be sufficient for heat removal, and since this 
size appeared reasonable for mobile applications, the amount 
of fuel was adjusted in order to give a radius within these 
dimensions. 

e. Moderator Proportion. The amount of moderator 
considered varied from zero to sixty per cent of the core 
volume. This had the effect of chenging the reactor from an 
essentially thermal one to an intermediate one, in which the 
greater part of the fissions take place at energies appre- 
clably above thermal. 

f. Reflector Composition. It was determined that 
as little as two per cent of the reflector volume was re~ 
quired to be sodium in order to carry off heat generated 
within the reflector. The amount of stainless steel was 
arbitrarily taken ag five per cent. 

By assigning the remainder of the reflector volume 


to the beryllium, it was found that the neutrons returned to 
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the core were largely thermal. In order to decrease the 
number of moderated neutrons reflected back to the core, ead~ 
ditional sodium was introduced into the reflector and the 
beryllium decreased, g0 that the percentage of beryllium 
considered ranged from a low of 55% to a high (one run only) 
of 95%. 

II. Effect of Various Parameters. 

ae Fuel Concentratione Apart from its effect on 
the critical radius of a reactor, the fuel concentration 
plays an important part in determining the energy spectrum of 
neutrons in the reactor. Neutrons are born fast and, during 
the slowing down process, may be absorbed by reactor ma~ 
terialse Absorption during slowing down increases the number 
which must be born fast in order that enough remain after 
Slowing down to cause further fissions and perpetuate the 
chain reaction. 

An examination of the cross sections of the ma~ 
terials considered here shows that the predominant absorber 
of neutrons having energies greater than thermal is the 
uranium. Increasing the fuel thus increases the flux in 
niga? energy grounse 

Actually there are three things which may happen to 
a neutron while slowing down: i1t may be absorbed; it may be 
lost tnrough leakage; or it may undergo scattering collisions 
and be slowed down. The fuel concentration is so 3mall 


that it can have little effect on the slowing down powere 
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The fuel's important effect is on absorption tendencies and, 
ag a consequence, on the leakage, since an increase in fast 
flux with more fuel increases leakage. 

It is important to note that of the neutrons abe 
sorbed by the uranium during the slowing down process, 
roughly four-fifths cause fission of the absorbing nucleus. 
Thus more fuel has the added effect of changing an originally 
thermal reactor toward a reactor in which fissions at 
energies above thermal predominate. 

Figure 1 demonstrates the above effect of fuel con- 
centration on the ratio of fast to thermal flux. The ratio, 
calculated at the reactor center, is plotted versus fuel 
concentration for twenty different reactors all having at 
least 30% beryllium in the core. 

It should be pointed out that the ratio of fast to 
Slow flux has an important effect on the ability of the 
reflector to assist in levelling the power distribution. A 
greater number of fast neutrons leaking into the reflector 
increases the number of slowed neutrons which the reflector 
can return to the core to compensate for slow leakage, which 
or vl otherwise depress the power at the edge of the core. 

b. Moderator Concentratione As mentioned above, 
the ratio of the fast to the slow flux depends on leakage, 
absorption, and slowing down. Changing the moderator pro— 
portion in the core affects the first and last terms. For 


simplicity of discussion its effect on leakage may be 
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neslected, if the central regions of the reactor are con- 
sidered. 

In comparing the remaining two terms--absorption 
and slowing, down--it is important to note that the criterion 
for comparison is the ratio of elther to the sum of the two. 
Tee 1S, tie total probability of one or the other hapopning 
(excluding leaxage) must be unity. For this reason, de- 
creasing on@ has the effect of increasing the other. 

tne beryllium moderator is the predominant con- 
tributor to tHe slowing down term py virtue of its hicn 
value of E Any decrease in beryllium tnus tends to @n- 
hance the absorption in uranium during slowing, down. Thus 
tne ratio of fast to slow flux varies inversely with the 
beryllium concentration as shown in Figure 2 where @acgt/ By 
is plotted as a fuaction of per cent beryllium for reactors 


of various compositions. 


Comemeact of fuel and »mod@rator coneentrations 
On tower Distribution 

As discussed above, eitner more fuel or less 
beryllium increases the ratio of fast to slow flux and thus 
raises the ability of th: reflector to return slowed neutrons 
to the core and assist in levelling the power distribution. 
Fige 3 shows this errect for various beryllium concentrations 
in an otherwise identical reactor. as the beryllium decreases, 


the power near the core-reflector boundary rises. 
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The uniformity of power distribution, measured by 


the retio of maximum to avel@ee power, deyends or the point 





at which aren occurs, thie may be either at tne center, the 
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IA should 


edge, or, ideally, ootn simudteneously. That is, 
Dee@aeiads lew CO cari tre cenver power in order to incresse 


es mucin ash possigle witnouc rimuliteaneously raising a ° 


ae. 
ave AX 


As an inaication of this, the retio of P 1 Bhs 5 


/P 
MAX" aver 


decresses to 1.16(as beryllium decreases from 30 to 20 per 
cent), then rises to 1.62 as beryllium is further decreased. 
in the lstter case the edge nower is limitian. 

Fie, 4 snovs a similar effect on the power dice 
tribution for verious fuel concertrsations. Figs. 4a and 4b 
tro voete Tivua vVerietions fcr various fuel coneentratiom. 
Big. 4b snows that with rore fuel the greater fast flux 
Perults in @ Miyerificant number of neutronr being returned 
to the core before beiz completely modersved. hese 21d 
$m ievetlim, The power iiistribution, 

Also in Pig. 4d (9.5% fuel) the ler;e ratio of 
fast to therm.1 flux results in some 15% of fissions 
taking, vince in the two higher energy FroupA, as contrasted 


, | ee erat ean, 


wilh only 9% for the reactor of Fig 4 
d. Keplector Composition 
ry . Fe : a é p + A koa é aoe =~ ia 
aewet tect of veorying the nropPeakicns of the tnree 
Pemiec tor comnponcniee wae invertigated. in genet, for tne re= 


flector comporiAtions studied, beryliium maker the lergest 


centribntion to the clowing Gon term and “he smallest to 








ee 


the absorption. Thereforc, decreasing the reflector's 
beryllium concentration effectively increases the aebsorn- 
tion tendencies while decreasing the mod@ration. Reflectors 
with less beryllium thus return fewer slowed neutrons to 
the core. 

Fie. 5, Which shows the power distributions for 
three reactors with identical cores but with different 
proportions of beryllium in the reflector, demonstrates 
this effect. The rise in power occasioned by neutrons 
returned from the reflector is greatest for the reactor 
having the most beryllium in the rerlector, least for the 


reactor With least reflector beryllium. 


@. neflector Thickness 

Figs. 6, 7, and 8 show tne spatial distribution 
Of ilux fer reactors with retPector thicznesses of 10, 20, 
and 50 cm respectively. Since fast neutrons which leax into 
Caeser vector gre Subject to a greater cnance of slowing 
down in a thick retlector, the thermal flux rises more 
in such a reflector. In turn, this causes a power rise near 
the core=-reflector interface, since mores thermal neutrons 
are returned to tne core. 

This provides a means of increasing the tnermal 
flux in the outer core region where the power normally is 
less than at the reactor center. By a proper choice of 


reflector this increased themsal flux in the outer core 


region can result in a power level as high as the power 
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at the reactor center. This represents the closest possible 
approach to a uniform power distribution for a given homo- 

geneous reactor, since the power distribution sags between 

these two maximums. 

Fig. 9 shows the effect of varying the reflector 
thickness on the power distribution for a given core. A 10 
em reflector gives a power peak at tne reactor center. 
nomewhere betwecn 10 and 20 cm reflector thickness the 
optimum power distribution for this reactor occurs. For 20 
and 50 cm reflector thicknesses a power peak occurs at 
the surface of the spherical core. With increasing reflec- 
tor thickness a limiting value of Pmax/Pavg is 
approached asymptotically. A plot of the neutron flux in 
Bigs SO fer the 50 ecm rerlector shows thet te Tlieein ali 
energy groups approaches zero asymptotically to the axis. 
Since the current density J = -Dd(NV)/dr (as given in 
Chapter i), a nearly zero slope to the flux curve means 
that there are very few neutrons escaping from the outer 
reflector surface. This indicates that a further increase 
in reflector thicxness would have little effect on the 
power distribution, so that this is essentially an 
infinite reflector. 

The positive slope (at the core-reflector inter- 
face) to the thermal flux curve means that J is negative, 
indicating a net flow of neutrons from the reflector into 
the core. It is noted that with reflectors of 10 and 20 


cm thicxness there is a net leakage of group III neutrons 
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into the reflector. However, with a 50 cm thick edn a Meee 
there is a positive slope to the group III flux at the 
interface, indicating a net return of group III neutrons 
from the reflector to the core. This tends further to in- 
crease tne power rise in the outer core regions. 

The significance of fissions in group III can be 
further investigated. Figs. 4a,4b, 6, 7, and 8 show flux 
plots for different reactors. Histograms in these figures 
show that approximately 60% of the total fissions take 
place in energy group III. Thus, one would expect that the 
spatial power distribution would follow the core flux plot 
of group III quite closely. An examination of the figures 
mentioned shows tnat the power distribution does follow 
the group III flux very closely except in the outer core 
regions where the thermal flux begins to rise. At this 
point the power curve also exhibits an upward trend. This 
reflects the fact that the thermal neutron fission cross 
section is approximately three and a half times that of 
the average group III cross section, so that a small inc- 
rease in thermal flux has great significance. 

Another trend which ean be noted from tne figures 
cit@a 16 that the group III flux distribution sags con- 
tinually as distance from the center increases. Thus 
power which depends primarily on group III will sag also, 
even though the thermal flux pel dis relatively flat. If 
more moderator is added so that the reactor approaches 
more closely a thermal reactor, the sag in the thermal 


flux in tne outer regions will not be compensated for 
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by a return of reflector moderated neutrons, since, with 
relatively few fast neutrons present there will be little 
fast leakage input to the reflector. rig. 4a illustrates 


this trend. 


f. Effect of Varying Na/Fe Ratio 


In Fig. 11 is plotted power distributions for 
three reactors identical except for the sodiun and stainless 
steel content of the core. The curves are very similar 
excent for the outer core region. 

Because Fe has a higher absorption cross section 
it is to be expected that increasing the core Fe content 
Will tend to raise the fast flux somewhat, thereby 
increasing fast leakage and tne number of slowed neutrons 
returnec from the reflector. Tnis point is supported by 
the plot of run (11) which has the highest value of power 
at the core-r fiector interface, and the hignest Fe con- 
tent. 

An additional effect is the variation in difrusion 
coefficient with varying Na/Fe ratio in the core. Since Fe 
has « higher scattering cross section (by about a factor 
of three) than sodium, and a higher atomic density (by a 
Similar factor) a greater re content tends to decrease D 
and thereby decrease leakage. It is thought tnat this 
partially explains the slight difference between the plots 
Oi Wns (1) and (10) @s shew in Pig. 11. Run (1), wath 


less iron and more leakage experiences a slightly greater 
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rise near the core-reflector boundary. 
since the two effects mentioned are in opposi- 
tion the above discussion is ratner hypothetical. In 
general, varying the Na/Fe ratio has little effect on the 


ewerali power distribution. 


&- Optimum tower Distribution 

in rig. 10 are plotted two optimum power distri- 
butions obtainec. Shapes of the curves, even though quite 
different core compositions were considere., are practi- 
cally identical. This tends to indicate tnat a limiting 
value of power uniformity nas been met. In other words, the 
power sag, observed evidentally cannot be further lessened 
for the reactor type and size considered. 

The two reactors considered in ig. 10 have 
Critical radii of 42 cm and 60 cm. Sine@ reactor size 
aifects the proportion of the core which may be»affected 
by reflector moderated neutrons, it is thought that size 
Will play an important part in determining uniformity of 
power. A reactor with a smaller core will have a larger 
proportion of vo.ume which experiences a power rise due to 
reflector moderated neutrons. However it can be seé€n in 
Fig. 10 that decreasing the size from 60 to 4e em has 
essentially no effect on power distribution. 

The two reactors considered in Fig. 10 have 
different material compositions and differcnt critical 
radii. They are representative of the class of reactors 
investigated. It can be concluded, therefore, tnat in 


the range of size and composition investigated the optinun 





ry 


ee 
power curve has a characteristic shape which is that 
shown in Fig. 10. The ratio of Pmax/Pavge for this curve is 


approximately 1.16. 


h. Comparison of Four and Six Group kethods 

rig. 135 gives a comparison of calculations by 
four group and six group methods for the same reactor. 
in general the results are in good agreement, since the 
Critical radii differ by less than four per cent. The ratio 
or fast to thermal flux was 3.4 for four groups and 4.7 
for six groups. These values are of the same order of magni-e 
tude. 

Tne plot of spatial power distribution in the 
core shown in Fig. 13 gives different values of Prax/Pave 
for the two methods, although in general the curves have 
Similar characteristics. Better agreement in power dis- 
tribution between methods seemed to exist for reactors 
having other compositions, although other runs on identi- 
cal reactors were not made due to time limitations. Bince 
both methods use the age diffusion equation (with minor 
differences wnich are covered in the discussion of results) 
the results should be very Similar. "=xperimental results 
for Na, Be, Fe intermediate reactors were not available 


for @eenec. On thie accuracy of the results. 





oe 
REFERENCES 


C. Be. Weber and H. H. Hirsch: "International 
Gonference on the Peaceful Uses of Atomic Energy", 
Volume 9, Paper 561. 


J; KR. Johnson ard Cy Es. Curtis# “Interfiational 
Conference on the Peaceful Uses of Atomic Energy", 
Volume 9, Paper 559. 


Reactor Handbook (Engineering), AECD-3646, U. 3S. 
Atomic Energy Commission, May 1955, (293). 


S. Glasstone: "Principles of Nuclear Reactor Eng- 
ineering", D. Van Nostrand Co., New York, 1955 (646). 


“a bib\ 





~8),= 
DISCUSSION OF RESULTS 


Tne principal results of tnis study have been the 
determination of criticul radius, neutron flux and rission 
power distribution, and tne ef-ect upon these quantities of 
varying the material compositions of the core and reflector 
as well aS varying, the reflector thickness. 

Tne metiod of attack involved modifying tne are 
diffusion equation so as to permit solutions to be ob- 
taaned by use of Whirlwind I digital computer. One metnod 
was developed in Cnapter I which used a finite difference 
aporoxi.iation to tne age diffusion equation with neutrons 
Givided into four energy sroups. A second zetiod was devel- 
oped in Chapter II wnicn used a method of numerical inte- 
@ration witn neutrons divided into six energy groups. 

(me results of calculations based uvon identical 
reactors using the two different metrods were in substantial 
asreement. llowever, «inor differences in spatial power dis- 
trivations and critical radius were noted for solutions by 
tne two netnods. 

Tne exolanation for differences in the radius and 
power distribution given by the two methods is thought to 
be the result of differences in calculated average cross 
sections for the various energy groups. The enersy groups 
are so wide that it is necessary to obtain an average cross 
section and diffusion coefficient for eacn energy group. Tne 


process of averaging, could result in incompatibilities 
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between average material pronerties calculated for the 
two methods, particularly in regions where the cross sec- 
tion changes rapidly with energy, as is tne case for the 
scattering resonance of sodium near 2700 CV. 

Another major difference between the methods 
involved tne determination of tne elastic scattering 
input frou energy group i to the next lower energy 
sroup itl. In tne four-gsroup method a factor f was chosen 
sucn tnat diave r ae The soiution was not sensitive to 
tne value chosen for f as shown by successive runs using 
improved values of f. Chapter I (vage 4) discusses the 
manner in which f can be improved from run to run on a 
tria. and error basis. Tne use of the f factor gives good 
results for wide ener, y groups. Tne disadvantage of 
naving to estimate f can be overcome by considering a 
greater number of enersy groups where tne variation of 
the slowing down density, q, over the group can be consi- 
dered to be linear. In this case sufficient accuracy 


in 


. . :é out 
could be obtained by the assumption that aq = S(q, +435 
VE 


ie 
In the six group methoc, however, the assumption 

is that the elastic scattering out of an energy group is 

ecual to the average value of the slowing down density 


ous 


in an energy Group, = q, ° This is equivalent to 


‘dave 
asSuming that the f factor ecuals unity. These two 

different assumptions concerning the elastic scattering 
out of an energy group probably account principally for 


any differences between results obtained by the two 
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methods. Tnis difference can be illustrated by the fol- 


lowing figure: 
Ce Se &..0 wo Dovini GT ainty 


a tame 


vs 
: Letharcy 
Bec Bh Aves 
— —< OUT_ 44 avg 
a 440°" # 
oe ee, Gr, ¢ 1 
Sa we ~Olid line represents 
@ s o 
AA De woe. FRuegsontanuons 
an are a 
0 606 10.81 12.20 (thermal) 


Lethorgy, wu 
Tne solid line represents the true distribution of slowing 


down density. Tne slowing down density, q, increases for low 

lethargy values corresponding to the birth energies (10 to 

Oi Mev) of fission neutrons. For lethargies greater than 

the birtn lethargy, a decreases due re ansorption and 

leaxage. ror the relatively wide lethargy groups snown it 

can be seen tnat the assumption used with the six group method 
ave . ~out 


that qy ont is not nearly as true as the representation 


of the true continuous slowing down curve by tne assumption 
tnat Cee = i Gees With a value of f chos@n so as to 
approximate the true curve as closely as possible. 

The disadvantages of using a relatively small 
number of energy groups are as follows: (1) the true 
elastic scattering out of an energy group is difficult to 
estimate correctly; (2) determining average neutron cross 
sections for a particular lethargy group introduces errors, 
particularly wnen cross sections vary rapidly with energy. 
Considering a larger nucber of groups greatly reduces 
wnese ,roblems. 


The following advantages result from using a 
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smali number of energy groups: (1) it is possible to 
neglect elastic scattering across groups (for the raterials 
considered); (2) the energy spectrum of rission neutrons 
need not be known with great accuracy if group widths are 
suc that all fission neutrons can he considered to be born 
on Ome* or two energy groups. 

mie four group prorpram igs arranged to include the 
@frects of inelastic scattering. however, in this inves- 
tigation inelastic scattering is neglected, since it is 
considered to be of only minor importence for group Widths 
and materials considered. 

The plots of neutron flux and power distribution 
gave consistent results indicating qualitative trends as 
parameters were variec. Tne parameters considered were; 

(1) fuel concentration; (2) reflector composition; (3) Na/®8e 
ratio; (4) Na/Fe ratio; and (5) reflector thickness. An 
erfort was made to choose reactors which would have the 
naximum power occurring equally at the reactor edge and 
ceater to get the smallest value of Pyux/Payg (ideally 

equal to unity). The best value that could be obtained was _ 
emo fOr ran (Cc). Pigs. 11 shows plots for two optimum 

power distributions. Cojparison of tnese two curves with 
other power distributions near the optimum seems to indi- 


cate that the general shape of the curve cannot be changed 


very much for tne range of Darameters investigated. However, 


these results would not apply to reactors having a much larger or 
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smalier critical radius since the shape of the optimun 
power curve would very lil:ely be considerably changed with marked 
changes in reactor size. 
A ratio of 1.16 ene Prax/Fave was the best value 


obtained for the reactor size and type discussed. 


The two programs which were developed for use with 
tne Wwnirlwind I digital computer operated very satisfactorily, 
and it is hoped that other investigators will find these 
workable programs to be helpful. Both programs were character- 
ized by rapid convergence to the answer, averaging only 
about four or five iterations before the convergence criteria 
were met. 


—_ 


For further work on tne sudject or uniform power 


0) 


distribution, it is Suggested tiat multi-region reactors be 
investigated. The Whirlwind programs developed in Chapters [I 
ana II can be modified so as to apply to reactors having 
several resions of uniform material properties. The programs 
would lengthen considerably, however, and it probably would 
be necessary to use auxiliary Whirlwind storage. 

It would be of interest to compare tne results of 
this theoretical investigation with experimental results in 
order to determine how closely a finite difference or 
numerical integration approximation to the age diffusion 
equation will represent reality. 


AN extension of this study to investigate reactors 
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of a different range in size using other moderators 


and coolants would supplement the results given here. 
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It was the objective of this thesis to study 
tne practicability of achieving a uniform power distri- 
bution in a homogeneous nuclear reactor, through return 
of moderated neutrons from a reflector to the core to 
compensate for leakage. Multi-group methods of solution 
of neutron diffusion equations were developed. While it 
was impossible to compare the theoretical results obtained 
With experimantal data, because of lack of such data, the 
results appear reasonable in that they demonstrate trends 
which are in agreement with physical considerations. 

Studies were hased on a single basic reactor 
chosen because of suitability to high output mobile 
applications. 

It was found, for the reactor type and size studied 
that, wnile complete uniformity of power could not be 
achieved, nevertheless it was possible to approach 


uniformity quite closely. The best value obtained for 


Prax/ Pave was 1.16. 
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APPENDIX A 


at 
For a typical power reactor, Glasstone gives an 
estimated volumetric power output for the core region of 


600 een, - It will be of interest to determine the ap 
rroximate power potential of the reactor of Run 8, which 
gave optimum power distribution with the maximum power 
points equal at the core edge and center. The critical core 
radius wag 60.00 cm. with sodium comprising 70 per cent of 
the core volume. | 

Because of time limftations, detailed heat transfer 
Calculations are not made. However, an estimate of the 
power potential of the reactor can be made based upon the 
following calculations: 

Core radius = 60.00 cm. = 1.97 ft. 

Core volume = 9.05 x 10° em 

Total power output at an assumed 600 watts/em- 

= 543 MV 

It will be assumed that the sodium coolant speed is 
5 ft./sec. For a mean temperature of 800°F, the sodium has 
the following Vee ie 

< =a5e 1b/ft? 

© = 0.31 Btu/lb F 
1--S, Glasstone, "Principles of Nuclear Reactor Engineering," 

D. Van Nostrand Co., Inc., New York, 1956, page 636. 


2-~"Reactor Handbook: Engineering," U. S. Atomic Energy 
Commission, August, 1955, (Geneva), page 266. 





aa? 


For power applications, it is hypothesized that a 
eylindrical reactor of 1.97 ft. radius can be chosen which 
will have a similar power aiatribution and volume to the 
spherical reactor. For 543 MW power output, the mean tem 
perature rise, AT, of the coolant flowing axially through 


70% of the reactor cross section is: 


Btu 
_ = 543 MW x 3.413 x 10° hr Ma 


yr... ibe, ss Ok ee 
5 te. x 3600° P82 x 52 aed % O32 ue x 0.7 x (1.97) ft 


= 546°R 

Thus a power output of 543 MW can be obtained from 
this reactor if sufficient heat transfer surface exists to 
give a mean temperature rise of Ch6 OF with a sodium coolant 
velocity of 5 ft/sec. With an inlet temperature of 600°F 
and a mean outlet eee of 1150°Fr, the maximum tem 
perature of a thin fuel element should not be excessivee 

The ability to transfer the required heat from the 
fuel elements to the coolant is the primary limiting factor 
which must be studied in determining the maximum power output. 

It is of interest to compare the output of this re- 
actor with the power requirements of a B52 tyne jet bomber. 
At 563 miles per hour, eight jets of 10,000 ibs. thrust each 
develop a total of 120,000 thrust horsepower. If an overall 
efficiency of 20% 19 assumed, this is equivalent to a heat 
source developing 450 MWe This can be compared to th 


roughly estimated reactor power output of 543 MW. 
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APPENDIX B 
TABS. I 


St 
MICROSCOPIC CROSS SzcTIONS 


A. S@attéring Cross Sections,o, (barns) 


Energy Lethargy, u= 
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O.L77% OIG D) WEG, S.5 4. Sey 
engl 4,606 |. 3.2 ESO 5.7 9.76 

pane. x 10°°S-1e5 ho7 4.6 1ox2 

coe OC 5.640 A. eet ©&,.0 107 

20.90 : 6.15 h,2 etn 10.9 

12.70 " 6.674 4.9 ee 5.9 10.9 
7.59 eT gee). re ws 8.8 Tages 
2oE5 “ 7.708 1%.5 Gee €.0 Jey 
2.69 4 G,225 200.0 ead 1639 
leo5 8.742 10.0 qa eG. O Lier 
0.9§5 "  G.259 4.5 9.3 Lame 
G.569 " 3.776 3.2 T.5 6.0 | Tee 
Gee ~" 10.293 oral eee) 20.6 
acne | 105010 8 ii ieee ~§ 64,0 19.2 

Pew 1Or° 12.142 Sel es 86G.O Ceo 

Vee ne Ak 4 3.1 ed 6.0 T6R0 
mo; 6 )~—CO 4. 806 oO Tet hO  10m0 
ig «= CCd16.138 3.2 ca oe “Seal 
OgeeS ~. 17.470 3.2 Jody: Se a eS 2.5 
QmO70 * 18.600 3.2 ites “Seg WV Ieno 


* ‘eutron Cross mmctopns, Biiles25, UmS. Atomae Energy 
Comission, Aug. 1955 (Geneva Suoplement) 
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PICRHOSCO; IC CROSS SECTIOVS (Continued) 


ry 


B. Absorption and *ission Cross Sections, o, andoy (barns) 

Gnergy Lethargy, u= U-2 
E (ev) 1n(10/E) Na re Pe of aa O 
10 0 O O O 1.29 1.55 O 
5.62 0.576 1.29 1.55 

7.16 1.151 O O O 1.29 1.55 O 
ares eT 1230 B2eGe 

© OO 2.303 O O O 1.27 1.52 O 
Oe562 2.879 L220 lee 

0.4316 % A544, O O O 1.30 1.56 O 
‘epalyaes 4.030 1.50 1.80 

GPL , 4606 O O O ie ome O 
6G.01lx 10-2 5.23 260° 2.0 
Bone) - 5.640 O O O 2.35 2.82 O 
20.90 6.157 ee00 E36 
Teer. 6.674 O O O 3.4 4,08 0 
1°00. Fmd OL 4.1 4,92 

i a at OO O 0 O moO 6.0 O 
B69 8,225 SOM 2 

i-cs. Serie 0 O O Pate | eo O 
6.955 : 9.259 oe “Ono 

569 : 9.776 O O O Lies toe O 

0.346 10.293 ty. hea 

O.-2CO0 =. 16,810 O O O ten “Oeee 9 
55.0 x 107° 12.142 52.0 62.4 

ie Mean O O O 55.0 66.0 0 
“07 -©% 18.806 200 46.0 

T7,O0y He 6.130 oh O.t ROO 662 4 O 
Cneco) s 17.470 0.15 0.1 O 50.0 187.5 
FOre7o.» Lescoo 0.36 C.6 C.GO06 320 375 0.0001 
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TAGE Li 
AVERAGE MLCRGSCOPIC CROSS SECTIONS ~ 


PER EWERGY GROUP 


A. Four Energy Groups (cross sections in barns) 


a a ee 


Group I (Lethargy range O to 4.C06) 


va Dial O is) 8 O O 
Be Set > 2.94 O O 
Fe 4.20 3.16 O O 
O Eeto 2.64 O O 
L=-5 6.64 6.64 12597 1.331 2.46 


Group II (Lethurgy range 4.606 to 10.81) 


Na Cees 2 lga6 QO O 
Be 5.94 eres 0) OQ O 
Fre Garo (eos QO OQ 
O Boas 5260 O O 
U=5 12.90 bee ygress 6.54 2246 
Group lII (Lethargy range 10.61 to 136.0) 
ua Soke j@0 ( 0.067 O 
pe 5.97 52 4 @) O 
re 1.4 ieeo O20C) O 
O Bye 3264 O O 
U5 a 0 L300 Gece CO.l 2.46 


0.0845 
0.209 
0.0359 
0.120 
0.0085 


0.0845 
0.209 
0.0359 
O«t20 
0.0085 


0.0045 
0.209 
0.0559 
Ogee ZO 


Group T (Lethursy 16.650, E = C.070 ev, 1000 F thermal energy) 





ila Bac Deed 0.30 O 
Be ore, EAC C.GO6 O 
Pe ities es HL hele C.60 O 
Geo 10.0 10.0 Bese a0. 2.46 


. Averaged with respect to lethargy. 
UeHd Stands for Uaas 


0.0545 
0.209 

0.0559 
0.0065 





-~lll- 


a 


CHOSS SECTIONS (Continued) 
Bb. OLX MMerey xroups 


6(fast): 0.234 Mev to 10 Mev 


5 : 5.5 kev to VU.23h Mev 

4 ; i256 CG to 5.5 kev 

3 : 3 ev to ioe ey 

2 : ~07 ev to 3 eV 

int pia APs MO? ey 
Na atomic density (x1U0<4)- eo 5k Dee fo!) 5 > , OMS 
Be « Var ey 1. > 209 
("@ - SCG > ~VYl119 . Se 
O(in Us) > .ULOL > .Unko : ee 
leew (90% ,in U0) oe 21, ©. W253 © wes 


Absorption Uross-v,cections (barns) (includes fission absorption) 


Energy Group Na Be Fe QO Ua235 Ym238 
ny ag) Uses omme.0 37i0u LG 
ba Cee Vane Use. 155.0 wee 
3 Oe Ome UO, UTE”) pis as 
L eo Ue Usws” UL Jape e eo 
S Ooo Ulte OU. wa. sme CAs 
6 Op UEw 0.0 .0 thee: es 
BiaStic Scattering Uross-3¢ctions (barns): 

1 303 lee 4 8 ee. yh 
iz Sire pry pile ies: re OSS 
3 Coe ae a oe aa5 ee) 
ib Pad re ac oe Te a ees et, 3 nO 
5 2.0 2 5e le wae 3. 5 grr, 8) 
6 Pes eee a ee ane aw OU 
rission Cross-sections (barns 

bb 32043 

4g Leos 

3 non 

ly LUE 

5 20 

6 cas: 





-l1l2- 
CROSS SECTIONS (Continued) 

Since BNL=325 lists total cross sections for 
elements other than fissionable materials certain 
assumptions are necessary in order to comoute scattering 
and absorption cross sections. The following assumptions 
were made with respect to all sree: other tnan 
U Hanan: 

(1) the thermal absorption cross sections listed 
in BNL-325 for 2200 m/s neutrons was decreased 
by a 1/v relation to the .07 ev thermal 
energy considered here; 

(2) the scattering cross section was assumed to 
be the total cross section at the point 
Where the absorption cross section, scaled 
by 1/v became negligible. 

Tnesé asfeunptions resulted in absorption in only the 
thermal and next higher groups. 

For the uranium the fission cross section is listed 
in BNL-3e5. The absorption cross section was synthesized 
from this under the assumption that the capture to fission 
ratio is constant. The scattering cross section of U-245, 
and all cross sections of VU-238, proved to be unimportant 
for this probler:. The capture to fission ratio of U-2355 


was taken to be 0.154. 





mn - WW N 


Xe) ~~ OV 


10 
11 
Le 


14 
1s. 
aK6 
17 
18 
AL, 


Es lig eo 


TABLE Tk 


REACTOR COMPOSI VIONS 


Core 


Pe 


30 
AN, 
30 
2O 
20 
2O 
30 
20 
10 
50 
30 
165 
cO 
30 
30 
60 
60 

O 
30 


Volume 


£@ 

Gf 
9.6 
Ba> 
9.6 
9.6 


Oat 
oa 
Oran 
oy 
29.7 
9.8 
9.6 
9.8 
9.8 
9.8 
9.8 
Oat 
9.7 


0.3 
O.4 
Os 
04 
0.4 
0.4 
0.3 
O.3 
0.3 
OE 
Oe 
Ome 
O.4 
O.2 
Ore 
O.2 
O.e 
0.3 
ORS) 


Fer Cent 
keflector 
Re 


20 
cO 
20 
10 
20 
40 
10 
10 
10 
20 


20 


40 
2O 
40 


UOs Na 


(Nn 


(es 


re 


OL 


ee 


a MH Ww WH OU WI WN 


Ur Ul Ul NPP SO 


OH 
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TABLE IV 
“eC EHLWIND INPUT DATA 
FOR x. OUR ENERGY GaouPs 


Energy 


Run ilo. Group D 4 FAs Za V2 
Cs) > I 1.867 0.756 0.02825 0.000106 0.000219 


ie 0.5686 1.412 O807646 0.000523 0.0Gia7 
pian 0.9769 1.44 0.05243 C.cCO8781l 0.01444 
Ay 0.9668 1.00 0.5196 0.03471 0.05247 


TR) I ie 0.68 0.06154 0.0 
Il Oe 251 gE 8) S: 0.12218 0.0 
el OMS 2D 1.05 O. 11559 0.000692 
Ey Greed 1.00 0.10581 0.00566 


2(C) if are oweal 0.50 0.02819 0.0001412 0.000291 
LT 0.569 eS. 0.07540 0.000697 0.00143 
EL 0.933 lee OROSe55 0.011136 0.01929 
bk 0.990 1.00 0.05186 0.04298 0.06909 


2(R) I 1.138 0.56 0.06154 0.0 
II 0.526 1.36 0.12236 0.0 
i =| (O.595 1.18 0.11539 0.000692 
70.601 1.00 0.10581 0.00460 


BiG) iE 1.370 0.50 O7.,0202/ 0.0U01768 0.060364 
stall Oesed eae 0.07540 0.040671 0.001790 

HE 0.982 great 0.05242 0.013490 0.02412 

af 0.986 1.00 0.05194 0.051228 ©,057e8 


oR) fe 1.138 O.56 ©.06154 0.0 
je 380.531 ies O.12296 0.06 

III 0.593 ies) Onllss9 0.000692 

7 0.601. 1.00 0.10581 0.00566 


4(C) i Cane Os6el O.cl473 O.0OCOlL412 0.00292 
188 Oe ee 1.528 0.065327 ©G.O00697 0.00142 
ital 1.189 Teee0 0.046105 0O.G11506 O0.O01gz 
ib 1.194 1.00 C.037016 0.043668 0.069¢9 


*(C) stands for reactor core properties, (R) for reflector 
Pepper ties. 





as 


WHIRLWIND INPUT DATA 
(continued) 


Run Energy 


No. Group _D__f £2, 2,  _VEe 
WR I 1.039 0.68 # £=0.06884 0 
7m ©.82 1308 093230. ) 


Laat 0.537 eae Me 0.12968 0.000522 
ra 0.544 1.00 0.12820 0.003914 
56 I rape AL 0.44 0.021473 0.0001412 0.000292 
tg 0.580 2.50 0.065327 0.000697 0.00143 
zit 1.189 1657 0.038105 0.011306 0.01929 
T 1.194% 1.00 0.037816 0.04367 0.06989 


5R ng 1.138 0.57 0.06154 0 
ia 0.531 2250 Oslez ie 0 

BE ony 0.593 200 0.11539 0.000692 

ay 0.601. 1.09 0.10581 0.00566 


bc  , —  -- -— | 


6C I 2.225 O«5 0.021473 0,0001412 0.000292 
It 0.580 2.0 0.065327 0.000697 0.00143 
Rit 1.789 1657 ©.038105 ©.012306 0.01929 
T 1.194% $1.00 0.037816 Q.043668 0.06989 
6R I 1.4004 0.85 0.04694 0 
TI 0.5445 1.22 #£0.10209 0 
TIE @.7501 1.02 O@0868 0.001033 


- o79e2 1.00 9.06587 0.006 
70 = I 1,867 0.621 0.02825 0.000106 0.900219 — 
Ii 0.568 1.528 0.07548 0.000523 0.0010 


Ligh 0.9797 1.660 0.05243 0.008781 0.0144 
7. 0.9868 paele 0.05196 0.03471 0.05247 
OR I 1.034 0.68 0.06944 0 
il 0.520 1.08 0.1331 0 
Tre = ate 1205 O51 3592 0.000491 
si 0.540 13060 0.1293 0.003773 
BC I 24209 0.621 0602158 0.000306 0.000219 


II 0.5796 1.528 0.96535 0.0005231 0,00107 
TIT 12.184 1.660 0.03813 0.005871 0.01444 


~  jeeo 1,00 0.03784 0.0354 0.05247 
BR = 1.039 0.68 0.06884 0 
II 6.522 1.08 0.13232 @) 


Lie O57 1.05 0.12968 0.000522 
T 0.544 1.00 0.12820 0.003914 


be a ee i | i | 
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WHIRLWIND INPUT DATA 











(continued) 
Run kbnergsy 
T \ 
‘No. Group D f oa 7 = vee 
9C I 2 700 0.621 0.01436 0.000f96 0.000219 
etd 0.591 1.528 0.05527 0.000523 6.00107 
Tl 8 1,498 1.668 0.02386 0.009121 0.01444 
T. 3660 1-00 0.02373 0.93631 0.05247 
9R it 1.934 0.68 0.06944 © 
II 02520 1.08 0.1331 ) 
III Gare 860.1 5 Oe 0.000491 
T Oe540 1,00 0.1293 0.003773 
106 I 1.698 0.9 0202919 0.00011 0.000219 
an 0.564 162 0.07287 0.000523 0.001073 
EY 0.779 1.4 0405522 0.00931 0.01444 
™ 0,724 ay 0.05474 0.03904 0.05247 
10R I leks DWC 0.06154 0 
II 0.531 095 0.12218 0 
III Ogos 860 0.11539 0.000692 
— 0,601 sas) 0.10581 0.00566 
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APPENDIX C 


Derivation of Equations 


Equation (5.5d) of Chap. I was stated without proof 
and used as a basis for the boundary conditions of the four 


group method. A proof of equation (3.5d) follows: 





Prove that 
7a re 
{F be Ny ~ —* 
— p+ c a - ee Fy 
s F ae eal NtOr MN ey ces 
—( Ytar Pyar) (365d) 


AS a preliminary step it is necessary to prove the 
following relationship 


7 + at +4 ar “ - 2 4, Nar 
Spode :o @ oo ars G,. 6 7a 
" ae (a) 


Proof 


Let = 4, +(s s — )( = 
(lam: ft Page + (EE) G 9 F 


k 
Cael 


a) 
aay Le 7 ay, 
7 








htt ae 
3 
_ Fn a = <a aa a t, a ae 
_ ae hf i és rn ~ Pr 3 
, hd/ h 
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APPENDIX C (continued) 


uSing tne relationship 


ce eae es 
nh h 


we have 


is r@dr = 2 f>[ 2, | y,,(2 oe hn) +, ra #25) 
_ P2 (4,034, (~or) te 14 (34, vor)] f 


Za era (a Lt Bat SG - 4 jar 


Thus, equation (a) is proved. Next introduce the variable Y 
and use the definition of VY in spherical coordinates 


5] 


HH 


o, 


r y 2 W mv i) 

; pM a ) 2 f je WE 
(~ wer sell = (or (ry eee - 
Ee } . 

Fe ii 
Use equation (b) in (a) with W replacing D , 
ER ne ae 


ec Y= (Dr )* »- Then 
oy W 
Wer Fp a SF (4, PHF Y,, 
( dp bh, Tat 2 hf ) 
Pry (yy, 27 + ¥ 4-1) - rn, (7 Ah + Yo) 
Dr Z i ae 


Vey = M0 I Fn 4, 24% t Yo 
y a 
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sek C (continued) 


N 
(-y+rs [5 ae n4i ra > APL, ae 


M NV 
“ay +2rey tard, ¥,. | 


‘ Vora Ye te Ye 
oe (y -2U, + Yy,) = A et 











Zar \ wt! 
br (% -24, +4) 1 Cee ee) 
- ee ee 
2 Ar | 6 
| f 7 _ WY 7 Y 
-; RY ne +nG@ tar( 
Ar vey vee Ie WG FOG aa 
Wy 4-8 Oy tre 
bw _ tye wt "Ny 
+ a ep vf >= Hg ary, Wy) G 
(Yaesu), (a-2%em) 
Vary) 6 
Meth Ar= ',,, -—h the result becomes (continued) 
nly y: AY, rr 4 / 
Ga Par Ww tM 
_y (v = Ww + y i= Yd) ~ 4 Aft Vy ny 
V4) V-~/ | 7a 
— 2Db 6 
HY, 1 nV Me 5 aw tM 
Yo = “4 





Ee Py } 2 ee 
— PAP (Yu) — Yue Ou i, N4, 





fa 4 Pewee 2 
Sea Ar Ga te 





finally 


A 
iC aie bh | Ae 
(- w+ 7. | adage (Vy, a i [ ty 


| ZOp 
f 


~i(Y i ¢ ) ee 3 j fh 44 4 
6 ol or Zar a 
since YW, is to be identified with f= 9 . This 1s the 
proof of equation (%3.5d) of Chap. I which is concerned 
Only with the first part of the above expression relating 


to evaluation at the point rye 
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